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(54) CONTROL DEVICE OF ENGINE 

(57)Abstract: 

PROBLEM TO BE SOLVED: To realize the optimization of a 
fuel injection amount, an intake air amount and EGR amount and 
improve an operation feeling and reduce an exhaust gas 
emission by constituting so as to carry out a fuel injection 
control, an intake control and EGR control synthetically in 
response to the operation of a driver. 

SOLUTION: In a control part 30, a target engine torque is set 
from an engine speed Ne and an accelerator opening S at 31 
and a basic fuel injection amount, EGR ratio and an equivalence 
ratio in a cylinder in response to this target value are set 
initially at 32. The equivalence ratio of a real EGR gas is 
estimated from the initial set value of the equivalence ratio in 
the cylinder and based on the estimation value of this 
equivalence ratio and respective initial set values of the 
equivalence ratio in the cylinder, the basic fuel injection amount 
and EGR ratio, the control target value for the partial pressure 
of the air component and non-air component of an intake pipe 
pressure are set at 34. The partial pressure of the non-air 
component is estimated by using an intake system model and 
EGR gas flow rate is set based on the deviation for the control target value of the partial pressure of 
the non-air component at 36. 
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CLAIMS 



[Claim(s)] 

[Claim 1] In the control unit of the engine which carries out adjustable control of fuel oil consumption and 
the throttle opening according to actuation of an operator A means to set up the desired value of the 
parameter according to an engine output state from an engine speed and accelerator opening, A means to 
initialize basic fuel oil consumption, an EGR rate, and the equivalent ratio in a cylinder based on the above- 
mentioned desired value, respectively, The equivalent ratio of actual EGR gas is presumed from the initial 
value of the above-mentioned equivalent ratio in a cylinder. The estimate of this equivalent ratio, The initial 
value of the above-mentioned equivalent ratio in a cylinder, and the initial value of the above-mentioned 
basic fuel oil consumption, A control-objectives value [ as opposed to the air component partial pressure of 
the pressure-of-induction-pipe force based on the initial value of the above-mentioned EGR rate ], A means 
to set up the control-objectives value over the non-air component partial pressure of the pressure-of- 
induction-pipe force, respectively, The above-mentioned non-air component partial pressure is presumed 
using the inhalation-of-air system model in consideration of the air component and non-air component in the 
air flow rate which passes a throttle valve, and EGR gas. The estimate of this non-air component partial 
pressure, A means to set up an EGR gas flow rate based on deflection with the control-objectives value of 
the above-mentioned non-air component partial pressure, The above-mentioned air component partial 
pressure is presumed using the above-mentioned inhalation-of-air system model. The estimate of this air 
component partial pressure, A means to set up the air flow rate which passes the above-mentioned throttle 
valve based on deflection with the control-objectives value of the above-mentioned air component partial 
pressure, and the air component in the above-mentioned EGR gas, A means to compute the control input of 
the actuator which adjusts the amount of EGR(s) based on the set point and the pressure-of-induction-pipe 
force of the above-mentioned EGR gas flow rate, A means to compute the control input of the actuator 
which adjusts throttle opening based on the set point and the pressure-of-induction-pipe force of an air flow 
rate of passing the above-mentioned throttle valve, The control unit of the engine characterized by having a 
means to set up the final basic fuel oil consumption for computing the control input to the injector which 
injects a fuel. 

[Claim 2] The control unit of the engine according to claim 1 characterized by calculating the estimate of the 
equivalent ratio of the above-mentioned EGR gas based on the output of an air-fuel ratio sensor. 
[Claim 3] The control unit of the engine according to claim 1 characterized by calculating the estimate of the 
equivalent ratio of the above-mentioned EGR gas as first-order lag over the initial value of the above- 
mentioned equivalent ratio in a cylinder. 

[Claim 4] The control unit of the engine according to claim 1 characterized by calculating the estimate of the 
equivalent ratio of the above-mentioned EGR gas by the first-order lag of the initial value of the above- 
mentioned equivalent ratio in a cylinder, and the pressure-of-induction-pipe force and the dead time set up 
by the engine speed. 

[Claim 5] The control unit of the engine according to claim 1 characterized by making into the measurement 
value by the inhalation air content sensor the air flow rate which passes the throttle valve in the above- 
mentioned inhalation-of-air system model. 

[Claim 6] The control unit of the engine according to claim 1 to 5 characterized by applying the value which 
carried out time quadrature of the error of the pressure response forecast of the above-mentioned non-air 
component partial pressure, and the estimate of the above-mentioned non-air component partial pressure in 
case the above-mentioned EGR gas flow rate is set up. 

[Claim 7] The control unit of the engine according to claim 1 to 6 characterized by applying the value which 
carried out time quadrature of the error of the pressure response forecast of the above-mentioned air 
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component partial pressure, and the estimate of the above-mentioned air component partial pressure in case 
the air flow rate which passes the above-mentioned throttle valve is set up. 

[Claim 8] The control unit of the engine according to claim 1 to 7 characterized by restricting to below the 
realizable maximum EGR gas flow rate in case the above-mentioned EGR gas flow rate is set up. 
[Claim 9] The control unit of the engine according to claim 8 characterized by setting up the above- 
mentioned maximum EGR gas flow rate according to a service condition. 

[Claim 10] The control unit of the engine according to claim 8 characterized by making the above- 
mentioned maximum EGR gas flow rate into a controllable value in 1 control period. 
[Claim 1 1 ] The control unit of the engine according to claim 1 to 1 0 characterized by restricting to below 
the realizable maximum air flow rate in case the air flow rate which passes the above-mentioned throttle 
valve is set up. 

[Claim 12] The control unit of the engine according to claim 1 1 characterized by setting up the above- 
mentioned maximum air flow rate according to a service condition. 

[Claim 13] The control unit of the engine according to claim 1 1 characterized by making the above- 
mentioned maximum air flow rate into a controllable value in 1 control period. 

[Claim 14] the above — the control unit of the engine according to claim 1 to 13 characterized by making 
final basic fuel oil consumption into the initial value of the above-mentioned basic fuel oil consumption. 
[Claim 15] the above — the control unit of the engine according to claim 1 to 13 characterized by setting up 
final basic fuel oil consumption based on the estimate of the above-mentioned air component partial 
pressure. 

[Claim 16] the above — the control unit of the engine according to claim 1 to 13 characterized by setting up 
final basic fuel oil consumption based on the pressure response forecast of the above-mentioned air 
component partial pressure. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the control unit of the engine which performs synthetically 

fuel-injection control, inhalation-of-air control, and EGR control. 

[0002] 

[Description of the Prior Art] Fuel quantity and an air content are determined recently as a basis of control 
of the engine output-shaft torque which is the physical quantity which carries out a direct action to control of 
a car. The technique of improving the responsibility over an operator's demand output and obtaining good 
performance-traverse ability is proposed variously. To JP,1-313636,A The technique which both carries out 
electronics control of the throttle-valve opening to if the desired value of engine output-shaft torque is set up 
from an accelerator control input and an engine speed and fuel oil consumption is controlled according to 
the set-up target torque with an actuator, and controls an inhalation air content is indicated. 
[0003] However, under the effect of the response delay of the inhalation air for being filled up with volume, 
such as an inhalation-of-air collector chamber of a throttle- valve lower stream of a river, at the time of 
transient operation, the effect of the delay of the equipment to which an inhalation air content is changed, 
etc. Since fuel oil consumption was controlled corresponding to target torque to a gap being between the 
demand air content corresponding to the target torque in the time of controlling the throttle valve, and the air 
content actually inhaled by the cylinder, it was difficult to supply the neither more nor less between a fuel 
and air. 

[0004] Phase lag compensation equivalent to the response delay of the actuator which controls the delay and 
the inhalation air content of inhalation air by collector chamber restoration of inhalation air was carried out 
to the fuel oil consumption itself as opposed to the target output torque to cope with this and the excess and 
deficiency of a fuel are prevented at JP,3-185248,A to cope with. 
[0005] 

[Problem(s) to be Solved by the Invention] However, in a Prior art, since it is the open loop control to which 
feedback is not carried out to the actual inhalation air content produced as an actuation result of a throttle 
valve, it cannot respond to the flow rate change by the bulb contamination produced with the flow rate 
deflection by the variation between the individuals of a throttle valve, the blow-by gas of the flow rate 
change by change of the opening area by the temperature change of a throttle valve, and a throttle valve, etc. 

[0006] Moreover, although carrying out phase lag compensation of a fuel or ignition timing corresponding 
to the response delay of a filled up part to the collector chamber of EGR gas or an EGR valve which passes 
an EGR valve for reduction of exhaust gas emission in the case of the engine which uses EGR together is 
also considered Since the responsibility of a throttle valve differs from the responsibility of an EGR valve, it 
is difficult to attain under a transient the set point based on the optimum value at the time of steady 
operation or the target air- fuel ratio which can be found from this set point, and an EGR rate. 
[0007] This invention was made in view of the above-mentioned situation, according to actuation of an 
operator, it is performing synthetically fuel-injection control, inhalation-of-air control, and EGR control, 
optimization of fuel oil consumption, optimization of an inhalation air content, and optimization of the 
amount of EGR(s) are realized, and it aims at offering the control unit of the engine which can aim at 
improvement in an operation feeling, and reduction of exhaust gas emission. 
[0008] 

[Means for Solving the Problem] In the control unit of the engine with which invention according to claim 1 
carries out adjustable control of fuel oil consumption and the throttle opening according to actuation of an 
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operator A means to set up the desired value of the parameter according to an engine output state from an 
engine speed and accelerator opening, A means to initialize basic fuel oil consumption, an EGR rate, and the 
equivalent ratio in a cylinder based on the above-mentioned desired value, respectively, The equivalent ratio 
of actual EGR gas is presumed from the initial value of the above-mentioned equivalent ratio in a cylinder. 
The estimate of this equivalent ratio, The initial value of the above-mentioned equivalent ratio in a cylinder, 
and the initial value of the above-mentioned basic fuel oil consumption, A control-objectives value [ as 
opposed to the air component partial pressure of the pressure-of-induction-pipe force based on the initial 
value of the above-mentioned EGR rate ], A means to set up the control-objectives value over the non-air 
component partial pressure of the pressure-of-induction-pipe force, respectively, The above-mentioned non- 
air component partial pressure is presumed using the inhalation-of-air system model in consideration of the 
air component and non-air component in the air flow rate which passes a throttle valve, and EGR gas. The 
estimate of this non-air component partial pressure, A means to set up an EGR gas flow rate based on 
deflection with the control-objectives value of the above-mentioned non-air component partial pressure, The 
above-mentioned air component partial pressure is presumed using the above-mentioned inhalation-of-air 
system model. The estimate of this air component partial pressure, A means to set up the air flow rate which 
passes the above-mentioned throttle valve based on deflection with the control-objectives value of the 
above-mentioned air component partial pressure, and the air component in the above-mentioned EGR gas, A 
means to compute the control input of the actuator which adjusts the amount of EGR(s) based on the set 
point and the pressure-of-induction-pipe force of the above-mentioned EGR gas flow rate, A means to 
compute the control input of the actuator which adjusts throttle opening based on the set point and the 
pressure-of-induction-pipe force of an air flow rate of passing the above-mentioned throttle valve, It is 
characterized by having a means to set up the final basic fuel oil consumption for computing the control 
input to the injector which injects a fuel. 

[0009] Invention according to claim 2 is characterized by calculating the estimate of the equivalent ratio of 
the above-mentioned EGR gas based on the output of an air- fuel ratio sensor in invention according to claim 
1. 

[0010] Invention according to claim 3 is characterized by calculating the estimate of the equivalent ratio of 
the above-mentioned EGR gas as first-order lag over the initial value of the above-mentioned equivalent 
ratio in a cylinder in invention according to claim 1 . 

[001 1] Invention according to claim 4 is characterized by calculating the estimate of the equivalent ratio of 
the above-mentioned EGR gas by the first-order lag of the initial value of the above-mentioned equivalent 
ratio in a cylinder, and the pressure-of-induction-pipe force and the dead time set up by the engine speed in 
invention according to claim 1 . 

[0012] Invention according to claim 5 is characterized by making into the measurement value by the 
inhalation air content sensor the air flow rate which passes the throttle valve in the above-mentioned 
inhalation-of-air system model in invention according to claim 1 . 

[0013] In invention according to claim 1 to 5, in case invention according to claim 6 sets up the above- 
mentioned EGR gas flow rate, it is characterized by applying the value which carried out time quadrature of 
the error of the pressure response forecast of the above-mentioned non-air component partial pressure, and 
the estimate of the above-mentioned non-air component partial pressure. 

[0014] In invention according to claim 1 to 6, in case invention according to claim 7 sets up the air flow rate 
which passes the above-mentioned throttle valve, it is characterized by applying the value which carried out 
time quadrature of the error of the pressure response forecast of the above-mentioned air component partial 
pressure, and the estimate of the above-mentioned air component partial pressure. 
[001 5] In invention according to claim 1 to 7, in case invention according to claim 8 sets up the above- 
mentioned EGR gas flow rate, it is characterized by restricting to below the realizable maximum EGR gas 
flow rate. 

[0016] Invention according to claim 9 is characterized by setting up the above-mentioned maximum EGR 
gas flow rate according to a service condition in invention according to claim 8. 

[0017] Invention according to claim 10 is characterized by making the above-mentioned maximum EGR gas 
flow rate into a controllable value in 1 control period in invention according to claim 8. 
[0018] In invention according to claim 1 to 10, in case invention according to claim 1 1 sets up the air flow 
rate which passes the above-mentioned throttle valve, it is characterized by restricting to below the 
realizable maximum air flow rate. 

[0019] Invention according to claim 12 is characterized by setting up the above-mentioned maximum air 
flow rate according to a service condition in invention according to claim 1 1 . 
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[0020] Invention according to claim 1 3 is characterized by making the above-mentioned maximum air flow 
rate into a controllable value in 1 control period in invention according to claim 1 1 . 

[0021] invention according to claim 14 — invention according to claim 1 to 13 - setting — the above — it is 
characterized by making final basic fuel oil consumption into the initial value of the above-mentioned basic 
fuel oil consumption. 

[0022] invention according to claim 15 — invention according to claim 1 to 13 — setting — the above — it is 
characterized by setting up final basic fuel oil consumption based on the estimate of the above-mentioned air 
component partial pressure. 

[0023] invention according to claim 16 — invention according to claim 1 to 13 — setting ~ the above — it is 
characterized by setting up final basic fuel oil consumption based on the pressure response forecast of the 
above-mentioned air component partial pressure. 

[0024] Namely, in this invention, if the desired value of the parameter according to engine output states, 
such as an engine torque, is set up from an engine speed and the accelerator opening by actuation of an 
operator Based on this desired value, basic fuel oil consumption, an EGR rate, and the equivalent ratio in a 
cylinder are initialized, respectively, and the equivalent ratio of actual EGR gas is presumed from the initial 
value of the above-mentioned equivalent ratio in a cylinder. The estimate of this equivalent ratio, Based on 
the initial value of the above-mentioned equivalent ratio in a cylinder, the initial value of the above- 
mentioned basic fuel oil consumption, and the initial value of the above-mentioned EGR rate, the control- 
objectives value over the air component partial pressure of the pressure-of-induction-pipe force and the 
control-objectives value over the non-air component partial pressure of the pressure-of-induction-pipe force 
are set up, respectively. 

[0025] The above-mentioned non-air component partial pressure is presumed using the inhalation-of-air 
system model in consideration of the air component and non-air component in the air flow rate which passes 
a throttle valve, and EGR gas. Next, the estimate of this non-air component partial pressure, While setting 
up an EGR gas flow rate based on deflection with the control-objectives value of the above-mentioned non- 
air component partial pressure The above-mentioned air component partial pressure is presumed using the 
above-mentioned inhalation-of-air system model, and the air flow rate which passes the above-mentioned 
throttle valve is set up based on the deflection of the estimate of this air component partial pressure, and the 
control-objectives value of the above-mentioned air component partial pressure, and the air component in 
the above-mentioned EGR gas. 

[0026] And while computing the control input of the actuator which adjusts the amount of EGR(s) based on 
the set point and the pressure-of-induction-pipe force of the above-mentioned EGR gas flow rate, the final 
basic fuel oil consumption for computing the control input of the actuator which adjusts throttle opening 
based on the set point and the pressure-of-induction-pipe force of an air flow rate of passing the above- 
mentioned throttle valve, and computing the control input to the injector which injects a fuel is set up, and 
fuel-injection control, inhalation-of-air control, and EGR control are performed synthetically. 
[0027] In this case, the estimate of the equivalent ratio of the above-mentioned EGR gas may be calculated 
based on the output of an air- fuel ratio sensor, and may be calculated by the first-order lag of the initial 
value of the above-mentioned equivalent ratio in a cylinder as first-order lag to the initial value of the above- 
mentioned equivalent ratio in a cylinder, and the pressure-of-induction-pipe force and the dead time set up 
by the engine speed. Moreover, it is good also as a measurement value according the air flow rate which 
passes the throttle valve in the above-mentioned inhalation-of-air system model to an inhalation air content 
sensor. 

[0028] Moreover, it is desirable to apply the value which carried out time quadrature of the error of the 
pressure response forecast of the above-mentioned non-air component partial pressure and the estimate of 
the above-mentioned non-air component partial pressure, in case the above-mentioned EGR gas flow rate is 
set up, and in case the air flow rate which passes the above-mentioned throttle valve is set up, it is desirable 
to apply the value which carried out time quadrature of the error of the pressure response forecast of the 
above-mentioned air component partial pressure and the estimate of the above-mentioned air component 
partial pressure. 

[0029] Moreover, in the value to which it was desirable to restrict to below the realizable maximum EGR 
gas flow rate as for the set point of the above-mentioned EGR gas flow rate, and this maximum EGR gas 
flow rate was set according to the service condition, or 1 control period, a controllable value is employable. 
[0030] In the value to which similarly it was desirable to restrict to below the realizable maximum air flow 
rate as for the set point of the air flow rate which passes the above-mentioned throttle valve, and this 
maximum air flow rate was set according to the service condition, or 1 control period, a controllable value is 
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employable. 

[0031] moreover, the above — the initial value of the above-mentioned basic fuel oil consumption may be 
used for final basic fuel oil consumption as it is, and it may set it up based on the estimate of the above- 
mentioned air component partial pressure, or the pressure response forecast of the above-mentioned air 
component partial pressure. 
[0032] 

[Embodiment of the Invention] Hereafter, the gestalt of operation of this invention is explained with 
reference to a drawing, drawing 1 - drawing 8 — the 1st gestalt of operation of this invention — being 
involved — drawing 1 — the block diagram of a fuel, inhalation of air, and an EGR control section, and 
drawing 2 — for the flow chart of a fixed manipulation routine, and drawing 5 , the flow chart of a fuel, 
inhalation of air, and an EGR control manipulation routine and drawing 6 are [ the whole engine control- 
system block diagram and drawing 3 / the flow chart of initialization routine, and drawing 4 / the 
explanatory view of an inhalation-of-air system model and drawing 8 of the flow chart of crank angle 
interruption routine and drawing 7 ] the explanatory views of gas column distinction. 
[0033] Drawing 2 shows the engine control system which performs synthetically fuel-injection control, 
inhalation-of-air control, and EGR control, and while the various sensors for detecting an engine operation 
condition are connected to the Maine control unit 20 which consists of a microcomputer which calculates 
various controlled variables, the various actuators for engine control are connected. 

[0034] As sensors connected to the above-mentioned Maine control unit 20 For every predetermined crank 
angle, a pulse signal The gas column distinction sensor 3 which outputs the pulse signal for the gas column 
distinction generated between the crank angle sensor 2 to output and the pulse signal outputted from this 
crank angle sensor 2, the accelerator opening sensor 4 which outputs the voltage signal according to the 
amount of treading in of the accelerator pedal which is not illustrated, There are the pressure-of-induction- 
pipe force sensor 5 which outputs the voltage signal according to the pressure of inhalation of air, the inlet- 
pipe temperature sensor 6 which outputs the voltage signal according to the gas temperature within 
inhalation of air, and air-fuel ratio sensor which detects air-fuel ratio 7 grade. 

[0035] Moreover, there is ignition coil 1 1 grade formed successively as actuators connected to the above- 
mentioned Maine control unit 20 by the injector 10 of each gas column which injects a fuel, and the ignition 
plug 12 for every gas column, and EGR valve 14 for carrying out adjustable [ of the throttle actuator 13 and 
the amount of EGR(s) for carrying out adjustable / of the throttle opening ] further is connected. 
[0036] As a function which computes the various parameters with which the above-mentioned Maine 
control unit 20 processes the signal from each sensors, and an engine operation condition is expressed It has 
[ whenever / gas column distinction section 21 and crank angle ] the pulse recurrence-interval time amount 
calculation section 23, the engine-speed calculation section 24, the accelerator opening calculation section 
25, the manifold total pressure calculation section 26, the gas-temperature calculation section 27 of 
inhalation of air, and the air-fiiel ratio calculation section 28 whenever [ judgment section 22 and crank 
angle ]. Further It has the fuel, inhalation of air, and the EGR control section 30 used as the center of engine 
control, and has each function of the injection pulse period calculation section 40, the fuel-injection- timing 
setting section 41, the injection pulse generating section 42, the ignition timing setting section 43, and the 
ignition signal generator 44 as a function concerning a controlled-variable output. 

[0037] Namely, the input configuration of the output pulse signal (crank pulse) from the crank angle sensor 
2 and the output pulse signal (gas column distinction pulse) from the gas column distinction sensor 3 
performs gas column distinction in the gas column distinction section 2 1 , and a location is judged 
[ whenever / predetermined crank angle / of the specific gas column which carried out gas column 
distinction / whenever / crank angle / corresponding to the crank pulse by which a sequential input is carried 
out ] in the judgment section 22 whenever [ crank angle ] by making a location into a criteria crank location. 
Moreover, whenever [ crank angle ], by the pulse recurrence-interval time amount calculation section 23, the 
input spacing time amount of a crank pulse is clocked, the elapsed time of a between is computed whenever 
[ predetermined crank angle ], and an engine speed Ne is computed from the elapsed time of 1 80-degreeCA 
in the engine-speed calculation section 24. 

[0038] Moreover, based on the output voltage value of the accelerator opening sensor 4, the accelerator 
opening (the amount of accelerator treading in) S is computed in the accelerator opening calculation section 
25, and the pressure-of-induction-pipe force (the sum of the air component partial pressure within inhalation 
of air and a non-air component partial pressure; manifold total pressure is called hereafter) Pm is computed 
based on the output voltage value of the pressure-of-induction-pipe force sensor 5 in the manifold total 
pressure calculation section 26. Furthermore, based on the output voltage value of the inlet-pipe temperature 
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sensor 6, the gas temperature Tm of inhalation of air is computed in the gas-temperature calculation section 
27 of inhalation of air, and an air- fuel ratio lambda is computed based on the output voltage of the air- fuel 
ratio sensor 7 in the air- fuel ratio calculation section 28. 

[0039] On the other hand, a fuel, inhalation of air, and the EGR control section 30 in a detail As shown in 
drawing 1 , the target set torque section 31, the 1st load and combustion-control manager 32, a load and 
combustion-control manager 32of ** 2nd a, the inhalation-of-air system multiplier calculation section 33, 
the inhalation-of-air control manager 34, the control-factor calculation section 35, the F/B control section 
36, and the electronics control throttle (ETC) directions section 37 ~ and Consist of the EGR directions 
sections 38 and the target engine torque Tei is set up based on an engine speed Ne and the accelerator 
opening S in the target set torque section 3 1 . When the basic fuel oil consumption and the EGR set point 
(EGR rate) corresponding to the target engine torque Tei are initialized by the 1 st load and combustion- 
control manager 32, by the inhalation-of-air control manager 34 The pressure desired value within 
inhalation of air is divided and set to an air component partial pressure and a non-air component partial 
pressure from basic fuel oil consumption and the EGR set point. By the F/B control section 36 The throttle 
actuator indicated value as a control input to the throttle actuator 1 3 which is a controlled system according 
to the following inhalation-of-air system models, While computing the EGR valve indicated value as a 
control input to EGR valve 14, the final basic fuel oil consumption for computing the control input to an 
injector 10 by the 2nd load and combustion-control manager 32a is set up. 

[0040] The flow rate Qa of the new temper of passing throttle- valve lb infixed in inlet-pipe la of an engine 
1 as the inhalation-of-air system model adopted by this invention is shown in drawing 7 (throttle passage air 
flow rate), The EGR gas flow rate (EGR valve passage quantity of gas flow) Qe which passes EGR valve 14 
infixed in 1 d of exhaust air reflux tubing from exhaust pipe 1 c to inlet-pipe 1 a is supplied in inlet-pipe 1 a. 
By expecting the air content of the part which is the inhalation-of-air system model it is supposed that is 
flowed into the cylinder of an engine 1 , and is filled up with the inlet-pipe volume by throttle passage air 
flow Q a and the EGR valve passage quantity of gas flow Qe The target torque set up from the accelerator 
control input and the engine speed is transitionally realizable on time. 

[0041] The air component within inhalation of air from the sum of new mind of passing throttle- valve lb, 
and the air component in the EGR gas which passes EGR valve 14 A part for the air which flows into a 
cylinder is removed. Throttle passage air flow Q a, If a gaseous equation of state is applied using the gas 
constant Ra of cylinder inflow flow Q so of the air component the EGR valve passage flow rate Qea of the 
air component in EGR gas, and within inhalation of air, the inlet-pipe volume Vm, gas-temperature Tm of 
inhalation of air, and an air component Time amount variation dPmo/dt of the air component within 
inhalation of air can be expressed with the following (1) types. 

dPmo/dt= (Qa+Qea-Qso) and Ra-Tm/Vm — (1) [0042] Moreover, the non-air component within inhalation 
of air (non-air component in EGR gas) The non-air component which flows into a cylinder is removed from 
the non-air component in the EGR gas which passes EGR valve 14. Similarly time amount variation 
dPmee/dt of the non-air component partial pressure Pmee within inhalation of air The gas constant Re of the 
EGR valve passage flow rate Qee of the non-air component in EGR gas, the cylinder inflow flow Q see of 
the non-air component in EGR gas, and a non-air component can express by the following (2) formulas. 
[0043] 

dPmee/dt= (Qee-Qsee) and Re-Tm/Vm — (2) The EGR valve passage flow rate Qea of the air component in 
the EGR gas in the above-mentioned (1) formula and the EGR valve passage flow rate Qee of the non-air 
component in the EGR gas in the above-mentioned (2) formula can be expressed like following (3) and (4) 
types, respectively by applying the equivalent ratio phi of the EGR gas in EGR valve 14 inlet port to the 
EGR valve passage quantity of gas flow Qe. 

Qea= (1-phi), Qe — (3) Qee=phi-Qe — (4) [0044] Moreover, the cylinder inflow flow Q see of cylinder 
inflow flow Q so of the air component in the above-mentioned (1) formula and the non-air component in the 
above-mentioned (2) formula can be expressed with the following (5) and (6) types using the number L of 
gas columns of the stroke volume Vs per 1 cylinder, volumetric-efficiency etav, and an engine, respectively. 

Qso=((Pmo-Vs)/(Ra-Tm)) -etav-(Ne-L/120) - (5) Qsee=((Pmee-Vs)/(Re-Tm)) -etav- (Ne-L/120) - (6) 
[0045] Therefore, if it replaces by the multipliers a, ba, and be which show the part in a formula in the above 
(1) and (2) types by the following (7) - (9) formulas with the application of the above-mentioned (3) - (5) 
type and a matrix type describes the above (1) and (2) types As shown in the following (10) types, the time 
amount variation of the air component partial pressure Pmo and the time amount variation of the non-air 
component partial pressure Pmee can express the condition within inhalation of air by the equivalent ratio 
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phi of throttle passage air flow Q a, the EGR valve passage quantity of gas flow Qe 5 and EGR gas. 



=(Vs/Vm) -etav- (Ne-L/120) - (7) ba=Ra-Tm/Vm 



d 

dt 
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- a 
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ba (l-d>) • ba 



<E> • b e J L<5 e 



Q a 



(8) be=Re-Tm/Vm - (9) 



■(10) 



[0046] By using the above inhalation-of-air system model, it is based on the time amount variation of the air 
component partial pressure Pmo within inhalation of air, and the non-air component partial pressure Pmee. 
Throttle passage air flow Q a and the EGR valve passage quantity of gas flow Qe are computable. In the F/B 
control section 36 Feed back the deflection of the desired value of the non-air component partial pressure 
within inhalation of air, and the non-air component partial pressure estimate which is the calculated value of 
a non-air component partial pressure, and the EGR valve passage quantity of gas flow Qe is set up. 
Furthermore, the deflection of the desired value of the air component contained in this EGR valve passage 
quantity of gas flow Qe and the air component partial pressure within inhalation of air and the air 
component partial pressure estimate which is the calculated value of an air component partial pressure is fed 
back, and throttle passage air flow Q a is set up. 

[0047] And while setting up final basic fuel oil consumption by the 2nd load and combustion-control 
manager 32a and outputting to the injection pulse period calculation section 40, in the ETC directions 
section 37, the throttle actuator indicated value Sa is set up from the manifold total pressure Pm and the 
throttle passage air- flow-rate set point Qa, and it outputs to the throttle actuator 13, and further, by the EGR 
directions section 38, the EGR valve indicated value Se is set up from the manifold total pressure Pm and 
the EGR valve passage quantity-of- gas-flow set point Qe, and it outputs to EGR valve 14. In addition, the 
inhalation-of-air system multiplier calculation section 33 and the control-factor calculation section 35 
compute the multiplier of an inhalation-of-air system model, and the multiplier of feedback control, 
respectively. 

[0048] In the injection pulse period calculation section 40, injection pulse period Tout as a control input to 
an injector 10 is computed from the basic fuel oil consumption Gf set up by above-mentioned fuel, 
inhalation of air, and EGR control section 30, it sets according to the fuel injection timing Tinj set up in this 
injection pulse period Tout and the fuel-injection-timing setting section 41 by whenever [ crank angle / of 
the specification which defined the injection pulse generating timer beforehand in the injection pulse 
generating section 42 ], and an injection pulse is outputted to an injector 10 to predetermined timing. 
[0049] Moreover, ignition timing Tig is set up based on an engine speed Ne and the target engine torque 
Tei, an ignition pulse generating timer is set by whenever [ specific crank angle / which was beforehand 
defined by the ignition signal generator 44 according to this ignition timing Tig ], an ignition signal is 
outputted to an ignition coil 1 1 to predetermined timing, and an ignition plug 12 is made to discharge in the 
ignition timing setting section 43. 

[0050] Hereafter, the fuel, the inhalation of air, and EGR control processing performed with the above- 
mentioned Maine control unit 20 are explained according to the flow chart of drawing 3 R> 3 - drawing 6 . 
In addition, as for the subscript added to each parameter, i expresses [ a initial value and * / desired value 
and (-k) ] that it is a value in front of k control period (it is a value in front of 1 control period with a 
subscript (-1)). 

[0051] When the ignition switch which drawing 3 does not illustrate is turned on, a power source is supplied 
to the Maine control unit 20 and a system is reset, It is the initialization routine by which interruption 
activation is carried out, and first, if CPU is initialized at step S10, control data will be initialized at step 
S20. At step S30 Inhalation-of-air system constants, such as the inlet-pipe volume Vm, the stroke volume 
Vs per 1 cylinder, the engine number L of gas columns, the gas constant Ra of an air component, and the gas 
constant Re of a non-air component, are set up, and it escapes from a routine. 
[0052] And after system initialization, while the fixed manipulation routine shown in drawing 4 is 
performed for every (every [ for example, ] 10ms) fixed time amount, interruption activation of the routine 
of drawing 6 is carried out for every crank pulse input. 

[0053] In the fixed manipulation routine of drawin g 4 , first, as processing of the accelerator opening 
calculation section 25, the accelerator opening S is computed by carrying out A/D conversion of the output 
of the accelerator opening sensor 4, and the manifold total pressure Pm is computed as processing of the 
manifold total pressure calculation section 26 at step S60 by step S50 by carrying out A/D conversion of the 
output of the pressure-of-induction-pipe force sensor 5. Furthermore, the gas temperature Tm within 
inhalation of air is computed at step S70 by carrying out A/D conversion of the output of the inlet-pipe 
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temperature sensor 6 as processing of the gas-temperature calculation section 27 of inhalation of air. 
[0054] When an air-fuel ratio lambda is computed by progressing to step S80 and carrying out A/D 
conversion of the output of the air-fuel ratio sensor 7, as processing of the air-fuel ratio calculation section 
28 subsequently, at step S90 An engine speed Ne is computed from the elapsed time of 180-degreeCA 
computed by the crank angle interruption routine of drawing 6 mentioned later as processing of the engine- 
speed calculation section 24. At step SI 00 The fuel, inhalation of air, and EGR control manipulation routine 
of drawing 5 are performed as processing of a fuel, inhalation of air, and the EGR control section 30, and 
basic fuel-oil-consumption Gf*, the throttle actuator indicated value Sa, and the EGR valve indicated value 
Se are computed on the basis of the target engine torque Tei. 

[0055] It progresses to step SI 10. Then, as processing of the injection pulse period calculation section 40 
Various correction terms and a reactive component are added and basic fuel-oil-consumption Gf* computed 
at the above-mentioned step SI 00 is converted into injection pulse period Tout. As processing of the fuel- 
injection-timing setting section 41 When fuel injection timing Tinj is set up with reference to the map which 
uses an engine speed Ne and the target engine torque Tei as a grid, at step SI 20 Ignition timing Tig is set up 
with reference to the map which uses an engine speed Ne and the target engine torque Tei as a grid as 
processing of the ignition timing setting section 43, and it escapes from a routine. 

[0056] Next, the fuel, inhalation of air, and EGR control manipulation routine in the above-mentioned step 
S 1 00 are explained based on drawin g 5 . In this routine, with reference to the map which uses an engine 
speed Ne and accelerator opening S as a grid as processing of the target set torque section 3 1 , the target 
engine torque Tei is set up at step SI 50, and the inhalation-of-air system multiplier calculation section 33 is 
processed at step SI 60. In this inhalation-of-air system multiplier calculation processing, first, while setting 
up volumetric-efficiency etav based on an engine speed Ne and the manifold total pressure Pm Based on an 
engine speed Ne and manifold total pressure desired value Pm*i in front of 1 control period (-1), volumetric- 
efficiency etav* to pressure desired value is set up. By engine-speed Ne, gas-temperature Tm within 
inhalation of air, volumetric-efficiency etav, etav*, and the inhalation-of-air system constants Vm, Vs, L, 
Ra, and Re The inhalation-of-air system multipliers a, ba, and be by the above-mentioned (7) - (9) type and 
the inhalation-of-air system multipliers ca, ce, and d by the following (1 1) - (13) types, and d* are 
computed. 

ca=a/ba=(Vs/(Ra-Tm)) -etav-(Ne-L/120) - (1 1) ce=a/be=(Vs/(Re-Tm)) -etav-(Ne-L/120) - (12) d =(Vs/ 
(Ra-Tm)) -etav - (13) d*=(Vs/(Ra-Tm)) -etav* - (14) [0057] As processing of the 1st load and combustion- 
control manager 32, based on an engine speed Ne and the target engine torque Tei, the basic fuel-oil- 
consumption initial value Gfi, the EGR set point EGRS, and the equivalent ratio set point faii in a cylinder 
are set up by refer to the map, respectively, and processing by the inhalation-of-air control manager 34 is 
performed at step SI 80 in continuing step SI 70. 

[0058] In processing by the inhalation-of-air control manager 34, the equivalent ratio estimate fai which 
presumed the equivalent ratio of the EGR gas in EGR valve 14 inlet port from the equivalent ratio set point 
faii set up previously first is calculated. And air component partial pressure desired value initial value 
Pmo*i, non-air component partial pressure desired value initial value Pmee*i, and manifold total pressure 
desired value initial value Pm*i are computed by the following (15) - (17) types from the equivalent ratio 
estimate fai, the equivalent ratio set point faii, the basic fuel-oil-consumption initial value Gfi, the EGR set 
point EGRS, inhalation-of-air system multiplier d*, and the target air-fuel ratio ABFT. 
Pmo*i = (1/d*) and Gfi-ABFT/faii - (15) Pmee*i= ((fai-EGRS)/(l-fai-EGRS)) -(Re/Ra) and Pmo*i - (16) 
Pm*i =Pmo*i+Pmee*i — (17) [0059] Although the above-mentioned equivalent ratio estimate fai can 
acquire the precision which was most excellent in using the equivalent ratio calculation value computed 
from the actual air-fuel ratio lambda when the air- fuel ratio sensor 7 is a broader-based mold air- fuel ratio 
sensor As shown in the following (18) types, in consideration of the transportation lag time amount of the 
combustion gas which occurs for piping of EGR etc., the equivalent ratio estimate fai may be computed by 
the first-order lag of the equivalent ratio set point faii with a weighted average from the equivalent ratio set 
point faii in front of k control period (-k). 

fai= (1-q) and fai(-l)+q-faii (-k) -(18), however q: Weighted average multiplier [0060] Although it is good 
also as a constant which set up the weighted average multiplier q beforehand when calculating the 
equivalent ratio estimate fai by the weighted average by the above-mentioned (18) formula, strictly Since 
the transportation lag time amount of combustion gas changes with service conditions, so that first-order lag 
can be set up the optimal by the service condition The weighted average multiplier q is set up from the 
manifold total pressure Pm, and, as for the equivalent ratio set point faii in front of k control period (-k), it is 
desirable to consider as the value in front of k period equivalent to the dead time set up with an engine speed 
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Ne and the manifold total pressure Pm. 

[0061] In addition, in simple, as shown in the following (19) types, the equivalent ratio set point faii may be 
set up as equivalent ratio estimate fai as it is. 

fai=faii - (19) [0062] Then, it progresses to step S190 and the feedback multipliers fl, f2, hi, h2, gl, and g2 
shown by the following (20) - (25) formulas are computed as processing of the control-factor calculation 
section 35 with the inhalation-of-air system multipliers ba, be, ca, and ce and the equivalent ratio estimate 
fai. 

fl=(l/(ba-dt)) -n - (20) f2=(l/(fai-be-dt)) -n - (21) hl=ca - (22) h2=ce/fai - (23) gl=g -- (24) g2=g -(25), 
however dt: Control period n : Weighting factor (0< n<l) 
g : integral control multiplier (g>=0) 

[0063] Subsequently, according to the inhalation-of-air system model which progressed to step S200 and 
was mentioned above, the F/B control section 36 which computes the throttle passage air- flow-rate set point 
Qa and the EGR valve passage quantity-of-gas-flow set point Qe is processed. In this processing, first, in 
order to presume each time amount variation of an air component partial pressure and a non-air component, 
the air component partial pressure model value Pfo and the non-air component partial pressure model value 
Pfee according to an inhalation-of-air system model are computed, next the air component partial pressure 
estimate Pmo and the non-air component partial pressure estimate Pmee are computed as each of these 
partial pressure model values Pfo and calculated value which doubled Pfee with the manifold total pressure 
Pm which is the actual measurement of the pressure-of-induction-pipe force. 

[0064] And feed back the deflection of the non-air component partial pressure desired value initial value 
Pmee*i and the non-air component partial pressure estimate Pmee which were computed by processing by 
the above-mentioned inhalation-of-air control manager 34, and the EGR valve passage quantity of gas flow 
Qe is calculated. Furthermore, the deflection of the air component partial pressure desired value initial value 
Pmo*i and the air component partial pressure estimate Pmo which were computed by processing by the 
above-mentioned inhalation-of-air control manager 34 is similarly fed back using this EGR valve passage 
quantity of gas flow Qe, and it asks for throttle passage air flow Q a. 

[0065] Specifically each partial pressure model value Pfo and Pfee Throttle passage air flow Q [ in front of 
the inhalation-of-air system multipliers a, ba, and be, the equivalent ratio estimate fai, and 1 control period ] 
a (-1), It is computed by the following (26) and (27) types using the EGR valve passage quantity-of-gas- 
flow set point Qe (-1) in front of 1 control period, the air component partial pressure model value Pfo (-1) in 
front of 1 control period, and the non-air component partial pressure model value Pfee in front of 1 control 
period (-1). 

Pfo=(l-a-dt), Pfo(-l) +(ba-dt), andQa(-l)+ (ba-dt) -(1-fai)- Qe (-1) -- (26) Pfee=(l-a-dt) andPfee(-l)+ (be- 
dt) -fai-Qe (-1) — (27) [0066] Next, using each model value Pfo and Pfee which were computed by the 
above (26) and (27) formulas, and the manifold total pressure Pm measured by the pressure-of-induction- 
pipe force sensor 5, the air component partial pressure estimate Pmo is computed by the following (28) 
types, and the non-air component partial pressure estimate Pmee is further computed by the following (29) 
types from this air component partial pressure estimate Pmo and the manifold total pressure Pm. 
Pmo=Pfo+ (Pfo/(Pfo+Pfee)) -(Pm- (Pfo+Pfee)) - (28) Pmee=Pm-Pmo - (29) [0067] And according to the 
deflection of the desired value of a non-air component partial pressure, and estimate, the EGR valve passage 
quantity-of-gas-flow initial value Qei is computed by the following (30) types using non-air component 
partial pressure desired value initial value Pmee*i, the non-air component partial pressure estimate Pmee, 
and the feedback multipliers f2, h2, and g2. 

Qei=h2andPmee+f2- (Pmee*i-Pmee) - (30) [0068] For a certain reason, that it is not necessarily an 
implementation possible value also saturates the range of the following (31) types (range below or more 0 
maximum-stream-flow (Qe) max), it is made into a controllable (implementation is possible) flow rate, and, 
as for the EGR valve passage quantity-of-gas-flow initial value Qei computed by the above-mentioned (30) 
formula, makes this flow rate the final EGR valve passage quantity-of-gas-flow set point Qe. 
0 <=Qe<=(Qe) max — (31) — although it is good also as a constant which calculated beforehand the above- 
mentioned maximum EGR valve passage quantity-of-gas-flow (Qe) max by experiment etc. in this case, 
since a controllable EGR valve passage quantity of gas flow is dependent on the manifold total pressure Pm, 
it can realize exact F/B control by using the value set up by refer to the map etc. based on the manifold total 
pressure Pm. 

[0069] When controlling an EGR valve passage quantity of gas flow, furthermore, a controllable flow rate 
(it can be made to change) Since it is restricted by the manifold total pressure Pm and the EGR valve 
passage quantity of gas flow Qe (-1) in front of 1 control period, The maximum EGR valve passage 
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quantity-of-gas-flow variation (deltaQe) max is set up from the EGR valve indicated value Se in front of the 
manifold total pressure Pm and 1 control period (-1). By using maximum EGR valve passage quantity-of- 
gas-flow (Qe) max computed by the following (32) formulas with this maximum EGR valve passage 
quantity-of-gas-flow variation (deltaQe) max and the EGR valve passage quantity-of-gas-flow set point Qe 
(-1) in front of 1 control period More exact F/B control is realizable. 

(Qe) max=Qe(-l)+(deltaQe) max ~ (32) [0070] Then, according to the following (33) types, the throttle 
passage air- flow-rate initial value Qai is computed using the above-mentioned EGR valve passage quantity- 
of-gas-flow set point Qe with the air component partial pressure estimate Pmo, air component partial 
pressure desired value initial value Pmo*i, the equivalent ratio estimate fai, and the feedback multipliers fl, 
hi , and gl . And the computed throttle passage air-flow-rate initial value Qai is saturated in the range of the 
following (34) types (range below or more 0 maximum-stream-flow (Qa) max), and the throttle passage air- 
flow-rate set point Qa is defined. 

Qai=hl, Pmo+fl and (Pmo*i-Pmo)- (1-fai), Qe - (33) 0<=Qa<= (Qa) max - (34) [0071] Also in this case, 
like the case of the above-mentioned maximum EGR valve passage quantity-of-gas-flow (Qe) max, the 
above-mentioned maximum throttle passage air- flow-rate (Qa) max is good also as a constant set up 
beforehand, and the value set up by refer to the map etc. based on the manifold total pressure Pm in 
consideration of the controllable flow rate may be used for it. Furthermore, the maximum throttle passage 
air-flow-rate variation (deltaQa) max may be set up with the manifold total pressure Pm and the throttle 
actuator indicated value Sa (-1) in front of 1 control period, and maximum throttle passage air-flow-rate 
(Qa) max computed by the following (35) formulas with this maximum throttle passage air-flow-rate 
variation (deltaQa) max and the throttle passage air- flow-rate set point Qa (-1) in front of 1 control period 
may be used. 

(Qa) max=Qa(-l)+(deltaQa) max - (35) [0072] By the above, if processing of the F/B control section 36 in 
the above-mentioned step S200 ends, it will progress to step S210 next, and based on throttle passage air 
flow Q a and the manifold total pressure Pm which were computed at the above-mentioned step S200, the 
throttle actuator indicated value Sa is computed by refer to the map as processing of the ETC directions 
section 37. Furthermore, as processing of the EGR directions section 38, based on the EGR valve passage 
quantity-of-gas-flow set point Qe and the manifold total pressure Pm which were computed at the above- 
mentioned step S200, the EGR valve indicated value Se is computed by refer to the map, and it progresses 
to step S230 at step S220. 

[0073] At step S230, basic fuel-oil-consumption Gf* final as processing of the 2nd load and combustion- 
control manager 32a is set up, and it escapes from a routine. As shown in the following (36) types, the basic 
fuel-oil-consumption initial value Gfi can be used for it as it is, this final basic fuel-oil-consumption Gf* is 
good also as control of fuel priority, and the engine torque which avoided the control delay of the fuel 
quantity by the detection delay of air, and followed accelerator actuation of an operator can be generated, 
and it can improve the responsibility over an operator's demand output. 

Gf^Gfi ~ (36) [0074] Moreover, since delay may generate [ as opposed to / in fact / a control-objectives 
value ] the response of the pressure-of-induction-pipe force in this case, as it is show in the following (37) 
types using the air component partial pressure estimate Pmo, the equivalent ratio set point faii, the 
inhalation-of-air system multiplier d, and the target air-fuel ratio ABFT, according to the realistic pressure- 
of-induction-pipe force, basic fuel-oil-consumption GP final in D-JETRO may be compute, and the Air 
Fuel Ratio Control nature of a transient can be improve. 

Gf |e =d-Pmo-faii/ABFT — (37) [0075] As opposed to the above fixed manipulation routine in the crank angle 
interruption routine of drawing 6 First, according to the number of the gas column distinction pulses from 
the gas column distinction sensor 3 generated between the crank pulses from the crank angle sensor 2, a 
current gas column is distinguished as processing by the gas column distinction section 21 at step S300. 
Furthermore, processing distinguish a gas column after [ whose ] following the number of the crank pulses 
generated succeedingly is performed, and distinction processing is performed [ whenever / crank angle ] at 
step S3 10 whenever [ by the judgment section 22 / crank angle ]. 

[0076] As shown in drawing 8 , with this gestalt, a crank pulse is outputted from the crank angle sensor 2 for 
every BTDC97" of each gas column, 65 degree, and 10-degreeCA. From a gas column distinction sensor 
Between BTDC97" of a 3 cylinder, and BTDC10" of # 1 cylinder which is a front ignition gas column, # 
Three gas column distinction pulses, # One gas column distinction pulse is outputted between BTDC97" of 
two gas column distinction pulses, #1, and # 2 cylinder, and BTDC10" of a front ignition gas column 
between BTDC97" of a 4-cylinder, and BTDC10" of # 2 cylinder which is a front ignition gas column. 
[0077] Therefore, the crank pulse whenever the gas column distinction pulse was inputted, after counting 
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the pulse number and inputting three gas column distinction pulses The crank pulse of BTDC97" of a 3 
cylinder, and the crank pulse after two gas column distinction pulses were inputted # The crank pulse of 
BTDC97" of # 4-cylinder, The crank pulse after one gas column distinction pulse was inputted is a 
BTDC97" crank pulse of # 1 cylinder or # 2 cylinder, and if pre- gas column distinction is # 4-cylinder and 
# 1 cylinder and pre- gas column distinction are # 3 cylinders, it will be distinguished from # 2 cylinder. 
[0078] Moreover, the crank location of BTDC65" and BTDC10" is judged with the pulse number from the 
crank pulse of BTDC97", the crank pulse of BTDC97" of # 1 cylinder is set to 0 (criteria location), 
whenever a crank pulse is inputted, it carries out sequential count-up with 1, 2, 3, and — , and it distinguishes 
a crank location according to the counted value from a criteria location. 

[0079] At continuing step S320, whenever [ crank angle ] as processing of the pulse recurrence-interval time 
amount calculation section 23 the elapsed time to this crank pulse interruption generating from the last crank 
interruption generating — that is The elapsed time to this crank pulse input from the last crank pulse input is 
clocked. The elapsed time for 93 degrees whenever [ from MT93 and the crank pulse input of BTDC97" to 
the crank pulse input of BTDC65" / crank angle ] the elapsed time for 32 degrees [ whenever / crank angle / 
from the crank pulse input of BTDC10" to the crank pulse input of BTDC97" ] It stores in memory by 
setting elapsed time for 55 degrees to MT55 whenever [ from MT32 and the crank pulse input of BTDC65" 
to the crank pulse input of BTDC10" / crank angle ], The sum total of each elapsed time MT93, MT32, and 
MT55 is used for calculation of an engine speed Ne as elapsed time of 1 80-degreeCA. 

[0080] At step S330, processing of the fuel-injection-timing setting section 41 and the ignition timing setting 
section 43 is performed, and fuel injection timing and ignition timing are determined. That is, while 
converting into injection timing the fuel injection timing Tinj set up by the fixed manipulation routine from 
the specific crank angle defined beforehand, the ignition timing Tig similarly set up by the fixed 
manipulation routine is converted into ignition timing from the specific crank angle defined beforehand. 
[0081] At step S340, when it is interruption of whenever [ specific crank angle / which this crank angle 
interruption defined beforehand as processing of the injection pulse generating section 42 ], an injection 
pulse generating timer is set. Further and at step S3 50 As processing of the ignition signal generator 44, 
similarly, when it is interruption of whenever [ specific crank angle / which this crank angle interruption 
defined beforehand ], an ignition pulse generating timer is set and it escapes from a routine. Consequently, 
an injection pulse is outputted to an injector 10 from an injection pulse generating timer to the injection 
timing determined at the above-mentioned step S330, a fuel is injected, an ignition pulse is outputted to an 
ignition coil 1 1 from an ignition pulse generating timer by the ignition timing determined at the above- 
mentioned step S330, and ignition by the ignition plug 12 is performed. 

[0082] Inhalation-of-air control and EGR control are performed, the fuel oil consumption corresponding to a 
target engine torque by the above — receiving ~ an air component and a ratio, presuming an air component 
The inhalation air at the time of being filled up with the inlet-pipe volume of a throttle-valve lower stream of 
a river, and the response delay of EGR gas, The effect of responsibility which compensates the hard 
response delay of a throttle system and an EGR system, and is different by the throttle system and the EGR 
system is removable. By furthermore, the thing for which feedback to the actual inhalation air content 
produced as an actuation result of a throttle system is performed The flow rate change by change of the flow 
rate deflection by the variation between the individuals of a throttle valve, and the opening area by the 
temperature change of a throttle valve, The flow rate change by the bulb contamination produced with the 
blow-by gas of a throttle valve etc. can be avoided beforehand, control responsibility can be improved, and 
the flattery nature to accelerator actuation can be improved. 

[0083] That is, while being able to realize optimization of fuel oil consumption, optimization of an 
inhalation air content, and optimization of the amount of EGR(s) and being able to aim at improvement in 
an operation feeling, and reduction of exhaust gas emission according to actuation of an operator by 
performing synthetically fuel-injection control, inhalation-of-air control, and EGR control, the 
controllabilities in an extensive air- fuel ratio not only including a SUTOIKIO field but the Lean field can be 
improved, and free Air Fuel Ratio Control according to a service condition can be realized. 
[0084] With respect to the 2nd gestalt of operation of this invention, drawing 9 is the block diagram of a 
fuel, inhalation of air, and an EGR control section, and drawing 10 of drawing 9 and drawing 10 is the flow 
chart of a fuel, inhalation of air, and an EGR control manipulation routine. 

[0085] As this gestalt is shown in drawing 9 to the 1st above-mentioned gestalt, it sets to a fuel, inhalation of 
air, and the EGR control section 30. While changing the F/B control section 36 of the 1st gestalt, and the 
contents of processing of the 2nd load and combustion-control manager 32a a little and being referred to as 
F/B control-section 36A, and the 2nd load and combustion-control manager 32b, respectively 2nd 
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inhalation-of-air control manager 34a is added for the inhalation-of-air manager 34 of the 1st gestalt to this 
1st inhalation-of-air control manager 34 as 1st inhalation-of-air manager 34. 

[0086] Namely, with this gestalt, it sets to 2nd inhalation-of-air control manager 34a. Non-air component 
partial pressure forecast Pmee* which is the theoretical pressure response forecast of a non-air component 
partial pressure, and air component partial pressure forecast Pmo* which is the theoretical pressure response 
forecast of an air component partial pressure are computed. By F/B control-section 36A In case the EGR 
valve passage quantity-of-gas-flow initial value Qei is computed, the time quadrature value of the error of 
non-air component partial pressure forecast Pmee* and the non-air component partial pressure estimate 
Pmee is used. Moreover, in case the throttle passage air- flow-rate initial value Qai is computed, he is trying 
to use the time quadrature value of the error of air component partial pressure forecast Pmo* and the air 
component partial pressure estimate Pmo. Furthermore, final basic fuel-oil-consumption Gf* is computed by 
the 2nd load and combustion-control manager 32b using air component partial pressure forecast Pmo*. 
[0087] With this gestalt, to the fuel, inhalation of air, and EGR control manipulation routine of the 1st 
gestalt shown in drawing 5 , as shown in drawing 10 While changing the contents of the processing in step 
S200 into processing by F/B control-section 36A from processing by the F/B control section 36 The 
contents of the processing in step S230 are changed into processing by the 2nd load and combustion-control 
manager 32b from processing by the 2nd load and combustion-control manager 32a. Between the EGR 
valve indicated-value calculation processing by the EGR directions section 38 of step S220, and the final 
calculation processing of basic fuel oil consumption by the 2nd load and combustion- control manager 32b 
of step S230, step S225 which processes 2nd inhalation-of-air control manager 34a is inserted. 
[0088] Although the air component partial pressure estimate Pmo and the non-air component partial 
pressure estimate Pmee are computed and the EGR valve passage quantity-of-gas-flow initial value Qei and 
the throttle passage air-flow-rate initial value Qai are computed after that like the 1 st gestalt in processing by 
F/B control-section 36A of step S200 Under the present circumstances, to the 1st gestalt, as shown in the 
following (38) and (39) types, the time quadrature value Imee of the non-air component partial pressure 
error in front of 1 control period (-1) and the time quadrature value Imo (-1) of the air component partial 
pressure error in front of 1 control period are applied, respectively. 

Qei=h2, Pmee+fi2, and(Pmee*i-Pmee) + g2 and Imee (-1) - (38) Qai=hl, Pmo+fl, and (Pmo*i-Pmo) - (1- 
fai) -Qe+gl, Imo (-1) - (39) [0089] Like the 1st gestalt, the EGR valve passage quantity-of-gas-flow initial 
value Qei computed by the above (38) and (39) formulas and the throttle passage air-flow-rate initial value 
Qai saturate the range of the maximum stream flow from 0, and are taken as the EGR valve passage 
quantity-of-gas-flow set point Qe and the throttle passage air-flow-rate set point Qa. 

[0090] On the other hand by processing by 2nd inhalation-of-air control manager 34a in step S225 First, the 
air component partial pressure estimate Pmo, the throttle passage air-flow-rate set point Qa, With the EGR 
valve passage quantity-of-gas-flow set point Qe, the equivalent ratio estimate fai, the time quadrature value 
Imo of the air component partial pressure error in front of 1 control period (-1), and the feedback multipliers 
fl, hi, and gl Air component partial pressure target correction value Pmoh* which is the pressure desired 
value equivalent to the set-up throttle passage air flow rate is computed by the following (40) types. 
Pmoh*=(l/fl) - (Qa+ (1-fai), Qe+ (fl-hl) and Pmo-gl, Imo (-1)) - (40) [0091] Furthermore, non-air 
component partial pressure target correction value Pmeeh* which is the pressure desired value which is 
equivalent to the set-up EGR valve passage quantity of gas flow with the non-air component partial pressure 
estimate Pmee, the EGR valve passage quantity-of-gas-flow set point Qe, the time quadrature value Imee of 
the non-air component partial pressure error in front of 1 control period (-1), and the feedback multipliers f2, 
h2, and g2 is computed by the following (41) types. 

Pmeeh*=(l/f2) - (Qe+ (f2-h2), Pmee- g2, and Imee (-1)) - (41) [0092] Subsequently, air component partial 
pressure forecast Pmo* is computed by the following (42) types using air component partial pressure target 
correction value Pmoh*, air component partial pressure forecast Pmo* in front of 1 control period (-1), the 
feedback multiplier fl, and the inhalation-of-air system multiplier ba. 

Pmo*=(l-fl and ba-dt) andPmo*(-l)+ (fl and ba-dt) -Pmoh* - (42) [0093] Moreover, non-air component 
partial pressure forecast Pmee* is computed by the following (43) types using non-air component partial 
pressure target correction value Pmeeh*, non-air component partial pressure forecast Pmee* in front of 1 
control period (-1), the equivalent ratio estimate fai, the feedback multiplier £2, and the inhalation-of-air 
system multiplier be. 

Pmee*=(l-f2 and fai-be-dt) andPmee*(-l)+ (fl and fai-be-dt) -Pmeeh* - (43) [0094] And while computing 
the time quadrature value Imo of the error of the air component partial pressure forecast Pmo* and the air 
component partial pressure estimate Pmo which were computed by the above-mentioned (42) formula by the 
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following (44) types, the time quadrature value Imee of the error of the non-air component partial pressure 
forecast Pmee* and the non-air component partial pressure estimate Pmee which were computed by the 
above-mentioned (43) formula is computed by the following (45) types. 

Imo =Imo(-l)+ (Pmo*-Pmo), dt - (44) Imee=Imee(-l)+ (Pmee*-Pmee), dt - (45) [0095] In simple, air 
component partial pressure target correction value Pmoh* by the above-mentioned (40) formula and non-air 
component partial pressure target correction value Pmeeh* by the above-mentioned (41) formula can 
mitigate the count load of CPU greatly, although it is also possible to consider as air component partial 
pressure desired value initial value Pmo*i and non-air component partial pressure desired value initial value 
Pmee*i and control precision falls a little, respectively, as shown in the following (46) and (47) types. 
Pmoh* =Pmo*i - (46) Pmeeh*=Pmee*i ~ (47) [0096] In processing by F/B control-section 36A In this 
case, the EGR valve passage quantity-of-gas-flow initial value Qei, When saturating the throttle passage air- 
flow-rate initial value Qai in the range of the maximum stream flow from 0, respectively, Respond to the 
size relation between a initial value and the maximum stream flow, set / clear the saturation flags E and A, 
respectively, and by processing by 2nd inhalation-of-air control manager 34a It may be made to set up the 
time quadrature value Imo of an air component partial pressure error, and the time quadrature value Imee of 
a non-air component partial pressure error according to the value of each saturation flag, and a count load is 
mitigable, securing control precision to some extent. 

[0097] That is, when the EGR valve passage quantity-of-gas-flow initial value Qei and the EGR valve 
passage quantity-of-gas-flow set point Qe are equal, the saturation flag E is cleared, and when the EGR 
valve passage quantity-of-gas-flow initial value Qei differs from the EGR valve passage quantity-of-gas- 
flow set point Qe, the saturation flag E is set. Moreover, when the throttle passage air- flow-rate initial value 
Qai and the throttle passage air- flow-rate set point Qa are equal, the saturation flag A is cleared, and when 
the throttle passage air- flow-rate initial value Qai differs from the throttle passage air- flow-rate set point Qa, 
the saturation flag A is set. 

[0098] And when both the saturation flags A and E are cleared, the time quadrature value Imo of an air 
component partial pressure error is computed by the above-mentioned (44) formula, and when either of the 
saturation flags A and E is set, as shown in the following (48) types, let the time quadrature value Imo of an 
air component partial pressure error be a value in front of 1 control period. Moreover, when the saturation 
flag E is cleared, the time quadrature value Imee of a non-air component partial pressure error is computed 
by the above-mentioned (45) formula, and when the saturation flag E is set, as shown in the following (49) 
types, let the time quadrature value Imee of a non-air component partial pressure error be a value in front of 

I control period. 

Imo = Imo (-1) — (48) Imee=Imee (-1) (49) [0099] Moreover, in processing by the 2nd load and 
combustion-control manager 32b of step S230, final basic fuel-oil-consumption Gf* is computed from the 
basic fuel-oil-consumption initial value Gfi according to the following (50) types using air component 
partial pressure forecast Pmo* computed by 2nd inhalation-of-air control manager 34a, the equivalent ratio 
set point faii, the inhalation-of-air system multiplier d, and the target air-fuel ratio ABFT. 
Gf*=d-Pmo* and faii/ABFT — (50) [0100] Although it is also possible to adopt processing by the 2nd load 
and combustion-control manager 32a of the 1st gestalt about calculation processing of this final basic fuel- 
oil-consumption Gf*, without using air component partial pressure forecast Pmo* by the 2nd load and 
combustion-control manager 32b of this gestalt Like this gestalt by predicting theoretically the response 
value of the pressure-of-induction-pipe force over the current amount of control operation using air 
component partial pressure forecast Pmo* According to the hard delay of a throttle system or an EGR 
system of operation, or the delay of processing computation time, the delay produced in an actual inhalation- 
of-air system is avoided, the effect of pulsation etc. can be removed, the flattery nature to a transitional 
change of an air content can be improved, and highly precise control of air-fuel ratio priority can be realized. 

[0101] Moreover, with this gestalt, since the EGR valve passage quantity-of-gas-flow initial value Qei and 
the throttle passage air-flow-rate initial value Qai are computed using the integral values Imee and Imo of a 
control error in case feedback to the actual inhalation air content produced as an actuation result of a throttle 
system is performed to the 1st above-mentioned gestalt, the flattery nature of desired value to disturbance 
can improve, and F/B control precision can be improved. 

[0102] For drawing 1 1 - drawing 13 , with respect to the 3rd gestalt of operation of this invention, drawing 

II is [ the block diagram of a fuel inhalation of air, and an EGR control section and drawing 1 3 of the whole 
engine control-system block diagram and drawing 12 R> 2 ] the flow charts of a fixed manipulation routine. 
[0103] As this gestalt is shown in drawing 1 1 to the 1st above-mentioned gestalt or the 2nd above- 
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mentioned gestalt, while adopting the consistency sensor 9 which replaces with the pressure-of-induction- 
pipe force sensor 5, and detects the consistency within inhalation of air Have the inhalation air content 
sensor 8 which measures a throttle passage air flow rate, and it sets to the Maine control unit 20. The 
manifold total pressure calculation section 26 which computes the manifold total pressure Pm based on the 
output of the pressure-of-induction-pipe force sensor 5 It changes into manifold total pressure calculation 
section 26A which computes the manifold total pressure Pm based on the output of the consistency sensor 9, 
and the throttle passage air-flow-rate calculation section 29 which computes the throttle passage air-flow- 
rate measurement value Qave based on the output from the inhalation air content sensor 8 installed in the 
inlet pipe is added. 

[0104] Moreover, as shown in drawing 12 , in a fuel, inhalation of air, and the EGR control section 30, it has 
F/B control-section 36B which processes using the throttle passage air-flow-rate measurement value Qave. 
In addition, although drawing 12 shows the example which changed F/B control-section 36A into F/B 
control-section 36B to the 2nd gestalt containing 2nd inhalation-of-air control manager 34a, it is natural. [ of 
the F/B control section 36 being changed into F/B control-section 36B to the 1st gestalt which does not have 
2nd inhalation-of-air control manager 34a ] 

[0105] In the fixed manipulation routine performed for every fixed time amount with this gestalt as shown 
in drawing 13 Between step S70 which deletes step S60 which computes the manifold total pressure Pm by 
carrying out A/D conversion of the output of the pressure-of-induction-pipe force 5, and computes the gas 
temperature Tm of inhalation of air, and step S80 which computes an air- fuel ratio lambda Step S71 which 
computes the manifold total pressure Pm based on the output of the consistency sensor 9, and step S72 
which computes the throttle passage air-flow-rate measurement value Qave based on the output of the 
inhalation air content sensor 8 are inserted. 

[0106] With namely, the consistency rho of inhalation of air which carried out A/D conversion of the output 
of the consistency sensor 9, and asked for it as processing of manifold total pressure calculation section 26A 
in step S71 , and the gas constant Ra of air and the gas temperature Tm of inhalation of air for which it asked 
at step S70 By the following (51) types, the manifold total pressure Pm is computed, and as processing of 
the throttle passage air-flow-rate calculation section 29 in step S72, A/D conversion of the output of the 
inhalation air content sensor 8 is carried out, and the throttle passage air-flow-rate measurement value Qave 
is computed. 

Pm=rho-Ra-Tm — (51) [0107] And to calculation of the F/B control section 36 of the 1st gestalt, or the air 
component partial pressure model value Pfo by (26) types in F/B control-section 36A of the 2nd gestalt, by 
processing by F/B control-section 36B of this gestalt, as shown in the following (52) types, it replaces with 
throttle passage air flow Q a in front of 1 control period (- 1 ), and the throttle passage air-flow-rate 
measurement value Qave based on the output of the inhalation air content sensor 8 is used. 
Pfo=(l-a-dt), Pfo(-l) +(ba-dt), andQave+ (ba-dt) -(1-fai)- Qe (-1) - (52) [0108] Although other processings 
are the same as the 1st gestalt or the 2nd gestalt, since they measure an actual throttle passage air flow rate 
and perform feedback control, they can improve control precision more with this gestalt. 
[0109] 

[Effect of the Invention] As opposed to the fuel oil consumption corresponding to the desired value of the 
parameter corresponding to engine output states, such as a target engine torque, according to [ as explained 
above ] this invention In order to perform feedback to the actual inhalation air content which performs 
inhalation-of-air control and EGR control, and is produced as an actuation result of a throttle system, 
presuming an air component and a non-air component, While compensating the hard response delay of the 
inhalation air at the time of being filled up with the inlet-pipe volume of a throttle-valve lower stream of a 
river and the response delay of EGR gas, a throttle system, and an EGR system The flow rate deflection 
remove the effect of responsibility which is different by the throttle system and the EGR system, and 
according to the variation between the individuals of a throttle valve, The flow rate change by the bulb 
contamination produced with the blow-by gas of the flow rate change by change of the opening area by the 
temperature change of a throttle valve and a throttle valve etc. is beforehand avoidable. That is, according to 
actuation of an operator, fuel-injection control, inhalation-of-air control, and EGR control can be performed 
synthetically, optimization of fuel oil consumption, optimization of an inhalation air content, and 
optimization of the amount of EGR(s) can be realized, control responsibility is improved, the flattery nature 
to accelerator actuation is improved, and the effectiveness excelled — improvement in an operation feeling 
and reduction of exhaust gas emission can be planned — is acquired. 
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WRITTEN AMENDMENT 



[ a procedure revision] 

[Filing Date] February 6, Heisei 9 
[Procedure amendment 1 ] 
[Document to be Amended] Specification 
[Item(s) to be Amended] 0048 
[Method of Amendment] Modification 
[Proposed Amendment] 

[0048] In the injection pulse period calculation section 40, injection pulse period Tout as a control input to 

an injector 10 is computed from basic fuel-oil-consumption Gf* set up by above-mentioned fuel, inhalation 

of air, and EGR control section 30, it sets according to the fuel injection timing Tinj set up in this injection 

pulse period Tout and the fuel-injection- timing setting section 41 by whenever [ crank angle / of the 

specification which defined the injection pulse generating timer beforehand in the injection pulse generating 

section 42 ], and an injection pulse is outputted to an injector 10 to predetermined timing. 

[Procedure amendment 2] 

[Document to be Amended] Specification 

[Item(s) to be Amended] 0058 

[Method of Amendment] Modification 

[Proposed Amendment] 

[0058] In processing by the inhalation-of-air control manager 34, the equivalent ratio estimate fai which 
presumed the equivalent ratio of the EGR gas in EGR valve 14 inlet port from the equivalent ratio set point 
faii set up previously first is calculated. And air component partial pressure desired value initial value 
Pmo*i, non-air component partial pressure desired value initial value Pmee*i, and manifold total pressure 
desired value initial value Pm*i are computed by the following (15) - (17) types from the equivalent ratio 
estimate fai, the equivalent ratio set point faii, the basic fuel-oil-consumption initial value Gfi, the EGR set 
point EGRS, inhalation-of-air system multiplier d*, and theoretical air fuel ratio ABFT. 
Pmo*i = (1/d*) and Gfi-ABFT/faii - (15) 
Pmee*i= 

((fai-EGRS)/(l-fai-EGRS)) -(Re/Ra) and Pmo*i 
-(16) 

Pm*i =Pmo*i+Pmee*i - (17) 
[Procedure amendment 3] 
[Document to be Amended] Specification 
[Item(s) to be Amended] 0074 
[Method of Amendment] Modification 
[Proposed Amendment] 

[0074] Moreover, since delay may generate [ as opposed to / in fact / a control-objectives value ] the 
response of the pressure-of-induction-pipe force in this case, as it is shown in the following (37) types using 
the air component partial pressure estimate Pmo, the equivalent ratio set point faii, the inhalation-of-air 
system multiplier d, and theoretical air fuel ratio ABFT, according to the realistic pressure-of-induction-pipe 
force, basic fuel-oil-consumption Gf* final in D-JETRO may be computed, and the Air Fuel Ratio Control 
nature of a transient can be improved. 
Gf*=d-Pmo-faii/ABFT - (37) 
[Procedure amendment 4] 
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[Document to be Amended] Specification 
[Item(s) to be Amended] 0082 
[Method of Amendment] Modification 
[Proposed Amendment] 

[0082] By the above, inhalation-of-air control and EGR control are performed to the fuel oil consumption 
corresponding to a target engine torque, presuming an air component and a non-air component. The 
inhalation air at the time of being filled up with the inlet-pipe volume of a throttle- valve lower stream of a 
river, and the response delay of EGR gas, The effect of responsibility which compensates the hard response 
delay of a throttle system and an EGR system, and is different by the throttle system and the EGR system is 
removable. By furthermore, the thing for which feedback to the actual inhalation air content produced as an 
actuation result of a throttle system is performed The flow rate change by change of the flow rate deflection 
by the variation between the individuals of a throttle valve, and the opening area by the temperature change 
of a throttle valve, The flow rate change by the bulb contamination produced with the blow-by gas of a 
throttle valve etc. can be avoided beforehand, control responsibility can be improved, and the flattery nature 
to accelerator actuation can be improved. 
[Procedure amendment 5] 
[Document to be Amended] Specification 
[Item(s) to be Amended] 0099 
[Method of Amendment] Modification 
[Proposed Amendment] 

[0099] Moreover, in processing by the 2nd load and combustion-control manager 32b of step S230, final 
basic fuel-oil-consumption Gf* is computed from the basic fuel-oil-consumption initial value Gfi according 
to the following (50) types using air component partial pressure forecast Pmo* computed by 2nd inhalation- 
of-air control manager 34a, the equivalent ratio set point faii, the inhalation-of-air system multiplier d, and 
theoretical air fuel ratio ABFT. 
Gf*=d-Pmo* and faii/ABFT - (50) 



[Translation done.] 
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1 aftfcl&iteStl. x>S» 1©->>;>^C^ 

QMiiMMQa. <t EGR^^U^jiii^fefiQ e i(Cj;-p 
0. T^-fe-iU«fpgix>^>|BiKSt*iP > ^i y /ceS 

h t> z mmmicmtitj: < mm? z c t w-c # 
[oo4i] ©mifrto^MfiS^-a. ^ h)i>'W7 

1 b^rilii-rsffmi. EGRAjl/^l 4*ilj§-rSE 
50 GR^t|i©^^S^<?:©fn3!P6. '> U > ^F*9^A-T 
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Qa, EGR^^O^a^ODEGR^Vl/^iiia^S R«f flO2g^01$0Sftld Pm o/d t 

Qea, iRmmftOimfmWDt' V > ySEASHIQ s OTOCDSC^-r C t^ti. 

d Pmo/d t=(Qa+Qea-Qso)-Ra • Tm/Vm CD 
[0 04 2] Sfc, »atrt©*SM^ (EGR^ ^R^*©*^«ffi»<DEGR^;U^ii3ftStt»Qe e x 
*©*a«fl») tt, EGR^i^l 4 EG EGR^X*0*S«JSS»©^y>y«A«*Qs e 

^HP m e e (DBf Fl^tt d P m e e / d t EG*10 [0 04 3] 

dPmee/dt=(Qee-Qsee)«Re - Tm/Vm —(2) 



*Qe e = <f> • Qe •♦•(4) 

[0044] s/c, ±Ecoa:«:*5waa«jaE»©jxy 

->V> ^AStQ s e e tt. ^ti-en, 1 arSS/c 0 



±IB CUSSCCtoW & E GR //^*©aflWE»©E GR;< 
;U^ilM^*Q e a t ±fB(2)^tCfc^-SE GR^'X* 
<D}$%MS&ft<D E G R^Vl/ ^ISaiQ e e E G R 
>^ yaiS^/XSSSQ e^CEGR^^^ l 4 ADtefctf 

Qe a = (I -<*>) • Q e —(3) ★ 

Qso = C(Pmo • V s )/(R a • Tm)) • 7? v (N e • L/ 1 2 0)-(5) 
Qsee = ((Pmee • V s )/(R e • Tm)) • t?v- CN e • L/l 2 0) 

-(6) 

[0 04 5] S£oT, ±iB(l) > (2)3;tc±IBC3)— COS 
a, ba. befl^M, ±12(1), C2)SC* v h ij 9 * 



a =(Vs/Vm)-7jv(Ne 
ba=Ra • Tm/Vm 
be=R e • Tm/Vm 



eiEGR^©SlJt4>i«CA0, R»«rt©tt!»* 

L/l 2 0) - (7) 

♦.(8) 
-(9) 



d 


Pmo 




— a 


0 




Pmo 




















dt 


Pme e 




0 


— a_ 




Pme e 





ba (l-0>) • ba 



0 • be 



Qa 



.Qe 



•(10) 



[0 04 6] ^(Dmsm^to&m^zc tit**) * 

Qa £EGR^)l7mM*fxmMQe £*Sffi-r £C t 
F/B$fJ^a53 6"Ctt. «^«P9<D#^^R» 

[0047] ^ur, mzcon^ ■ jB»^ffl!-**-^* 
3 2 at«fB4S#«SBWiI*^Lr«i^^ 
B*IB»H«4 0^ffiarS±±fc^ ETCJg^SP3 7 

* h;U7**:*X-* 1 3^tH#U 36 50 



EGR^a$3 8-C, v^^;l/K±EPmiEGR 
^^U^iljft^SMIRSfflQ e Lfrh E GR^Jl/^fg^ 
ffiS eiSSltEGR^l 4^fcH;*rr£ 0 j^J, » 
AKffittfl[US3 3. «mfl%ffm3 5«, *ft* 

ft, awuR-t^KoflHK, ? - hv<* fM'&cDGm* 

[0048] *«*^<^^«FlBimi»4 0 rtt. ±BE«» 

• ®m • e g r $fj$BJ sis 3 o v&mutcm&iimmm&G 

;^«FffflTou tSrfltfflb, CCDPBW^U^^pfaTo u 

t &MRM(KS94 1 x^MLtcmmmmT i n j t 

A S>ifViWtt'WX&'t>Zsx#9 i o^m^-r&o 
[0049]^, J**BBBKJE»4 3 T?tt, x>J?> 

'A&fMT i gr*»5£U j£*{I-«H6£S|54 4r, 
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[0 0 5 0 ] ±.mJJ>ffl'&*~v bZOtCjz-o 

i/v ^t*^^s3o-c, ^atgsvm, i aerate oo 

[0 05 2] ^LT, ^rA^-i/t7>fX^ v 04 
0C^*rSJ8«!ffl;U-^>3&J— S^fHIS (^i* 1 Om 

[0 05 3 ]i4©S«l;^f>tlt ST. Xf 
2 6C[)MS<!:OT, »aeEZ3-fe>1f 5(Dffl*4A/D 

[0 054] Xf * t/SSO^jt^ SdStbfittH 

^2 8<7)#t3<hLT, 2mtk*>1>-7<Dtiijl*:A/D& 

*orffliisttA*inB'r4i, xx^^ssor, x>^ > 

ca=a/ba=(Vs/(Ra 
ce=a/be=(Vs/(Re 
d =(V s/(Ra • Tm)) 
d* = (V s/(Ra • Tm)) 
[0 057] ffi<Ar~v7si7CrC\Z > f 1©M 

WMTflBSWEiiGf \ E GRi£3!fiIE GRS, f>V> 

yrtSBJttsseffi f a i 1 ^ ^n-en-7->^*#M^ct 
9«5£U xf y^si8orsR5l*iJi3pv^-^i,3 4tcj: 

[0 05 8 ] KfU|lJ3Pv*--^ + 3 4CCcfc*fiifflr« % 

sr. 5fe«caafleuteaaitH3effl!f a i *^6egr^ 

Pmo* 1 =(l/d*) - Gf 1 • 
Pme e fc< = 

C(fai • EGRS)/(1 - f a i 
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*>@^#?ffi^2 4<D^a<tLT, fia-T4BI6<D^^ 
>^ftf«jA^;U-^>-CSttl$ntel 8 0" C A ©gift 
«FlH**6x>^>@e»Ne%»ao % x-r^^sioo 
~C\ l^-Ka-EGR«[IBlSB3 0©«HHiL/Tia5CD 
*R» • QRtl • EGRWawm^-^>%Hff 0, B«x 

f\ XD ? h Jl/T ^?^x-^}g7nffis a, EGRa* 

[0055] -ecofi, r/suo^jt^ pgw^^i/^ 
10 B$r^maigP4 0(DMa<bLr. iiexr ^^sloo-cjuu 

t£S3$4 1 CD^S<b It, x>^>@iSKN e <b g^x 

ESfJKT i n j *W0£?2>t. Ztv7'S12QV, **B$JtB 

>>> h-Jl/^T e ' £ i -T £ -7 * «#«UTfiife 

[0 05 6 ] ±IBX r * ^SLOOCCteW &#M4 ■ R 

20 EGRWiSP«ra^-^>tC^l»rBI5K:Sr5*Si58 

i8t§B3 1 <D%mt Utx>^>§®N e iT^tJl/ 

;u ^ T e 1 *«s u ^ ^sieor gRM^^^ffigp 

3 3<d#hs**t5. c©K^flHWHW«rw: % £ 

x>*P>IeH@&N eiv^W K£EEPm<LK:Sr3 

i«cl|rJCiTEE*a««K:»TS^«i*7jv *R3E 
30 U x>^>|5flBRNe, KSiertO^iaaTm, » 
a«fi*7?v f 7?v\ RSl^SVm.Vs.L.Ra.ReK: 
£0, f?MOC7)-(9)Si(Cj:^®^^a , ba.be, 
&2>\ «TO013--C13)a«:cfcS«»UR«BKca f ce f 

Tm)) - i?v(Ne -L/12 O)-(H) 
Tm)) • ?7 v (Ne -L/12 0 ) • (12) 



7?V 

7? v* 



-(13) 
-(14) 

Kfrzfl 4APCCteW4EGR^OS*tk*«jeLteS 
40 MttJiSffif a i*3fc»S. *lt. aStbJi^fflfa 

i . astbs^fii f a i \ m*imm*taxmtga£a 

Gf\ EGRtSfEGRS, OttURfltKd\ SUm 
«6ttABFT*6, «T©C15)^C17)SCCJ:0, 



ABFT/f a i 



(15) 



E GR S )) * (Re/Ra) - Prao' 
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Pm' 1 = Pmo* 1 + Pme e' 1 

[0059] ±^mtmmm faiit mmt^v- 

£f#SC<h*n?££>^ «T<D(i8)^t^TJ:Mc> E* 
f ai=(l-q) - fa i,. 15 + q 

x-mmicmmv mmw-^m q £ ^- *;u f 

f 1 = (l/(ba- dt))-n 
f 2 = (l/(f a i • be- d t )) 
h 1 = c a 
h 2= ce/f a i 
gl = g 

d t : mmmm 
mfrom ( o < n < 1 ) 
mftmmmtt ( g ^ o ) 

/BtaiaaiJ3 6©«iffl*tT^. cco^itu, ^ 

afr> B^^^r^tcSe-Dfc^^^EE^^l/ffiP f o 
RO'^^JS^E^e^HBP f e e^lTEHU 
cn^CD^EE^-r^ffiP f o . P f e e *5R«fflE27 
CDHiBfit?*57^*;l/ F^JEPm^C^#/cltSMi 

SlPmee^SffltSc 

[oo64]-eur v «rjso»«*iw^*-2>.|.34cc* 

P f o = (1 -a • d t) • P f oc- 
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-(16) 
"•(17) 

<t v % m&tmmm t a i 'o-rrawvcsat 
Jt^fflf a i iwmox^m^o 

f a i \_ k1 -(18) 

* [006 i]^ nawcctt. «TO(i9)S;cc^-rj:^ 

fc, 3JUtH5E«f a *©**S«t«lffi«f 
10 aiiltif ltfc^ 0 

f a i = f a i 1 — (19) 

[0 06 2] ^<D&, Xf'; ^S190^il*, ffllf&tigfcff 
ffl»3 5©«HliL/r, RMSMba.be.ca.cei 
S*it*^ffi f a i 4&CJ;9, JaT©(20)-(25)SerC?S 
t7-f-F^7^«Kf 1, f 2, hi. h2, gl, 

*2*»a-rs e 



fib 

n 



-(20) 

n .-.(21) 
-(22) 
"(23) 
♦(24) 
-(25) 

★Jr^MirgHJ L3fc*2«UJE^BEB«it»aK8SttP 
me e' 1 <b #S£ttS&&E iSSffi P m e e tcom^y 
j - F^ y^LXE GRs*)Vz7MM*fXifcMQ e 

CCDEGR^l/^ffij&tfx^SQe^rffi 

TjltfJ L /c^^^JEa^fif^^ti p m o ' 1 4 m 
30 ^S^JEJi^MPmoicDtB^^y F'<?*U 
XUv FJUiBiSSaiSiaiQ a 

[0 06 5 ] '=Rfl#JK:tt v S4WE*t^M8P fo, Pf 
eet*. R£UR«»a ,ba,b e, ^fiJtJt^ffi f a 
i . 1 PmmMmcox n y h ^ili^SiifcfiQ a c _ a , , 

0(-i,, 1 WKBaaSfflro^^SflJiR^EE* fvHB P f e 
e c -n», «T©C26) f C27)SC0Cj:D#fflSna Q 
n + (ba- d t) • Qa ( . 1} + 



Pfee = (l-a-dt) 



[0 0 6 6 ] ±iB(26) t (27)^r#aU/c^7 r 
JHBPf o. Pfeei v ®^©BE^*fe>1f 5-CltS'JU 
F±JEP mi*ffll\ «T<D(28)S;«:cfc «9 £A 
Pm o = P f o + 

(P f o/(P fo+Pfee)) 
Pme e = Pm — Pmo 

[oo67]-eu, #a^»»EEa«tt«nwa93effl so 



(ba-dt)-(l-fai)-Qe ( . n -(26) 
P f e e c . T) + 

(b e • d t ) • f a i • Qe c _ 13 -(27) 

^M^^ElfetffiP m o £|?£ti U . £6$c, 

*}ffKlPmo £-?~7fr;l> FiEPmi^6, KT 
<D(29)5£CC J: ^ifStKOt^smeUPm e e 



(Pm-(P fo + Pf e e))-(28) 

-(29) 

Pmee*\ ^^El^fiPmee, 
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A^ft»f2. h2, g2», tTF<D(30)5£CC<J: * E GR^';^I»^xSlR]l5S£iQ e 1 ^Sfflt 

Qe' = h2 • Pmee+f2 • (Pme e^-Pmee) -(30) 



[0 0 6 8] ±nEOO)^rCWmLtcEGR^)lzfmMtf 

&t>C£fe&£/t#>, «T<D(31)5£<D^H (OtLbS* 
^S(Qe)maxJ^TCDffiB) ifrM®JpJ#g (SI 

O^Qe ^(Qe)max -(31) 10 

/B^£^ir£C£^#£ 0 
[0 06 9] EGR^)l-7mj&rtXffi.m*ffl 

EGR^I/^ffiMtf^i&fiQec.,, tlc£^XfflmZti 20 
&/ctf>, -?~*)\sF±KPm± imWmMmCDEGR^X 

Qa^hl - Pmo+f 1 • (Pmo' 1 



M^{ta(AQe)max€r^L/, C(D^E GR^;^ 
Mii^^fi^{tS(AQ e )maxi 1 ffl'&MMm<D E G 
Rvs;l/^3lM^^d£fi^fiiQ e £ . a3 <b*c J: oTtTFO 
(32)5£T'l?ffi L ft:** E G R / W Si(Qe) 

(Q e )max= Qe MJ + (AQ e )max '(32) 

[0070] ^cDfg, jtiaE GR^t;i/^aa^ij«» 

#IEg^fiI*7J»il££tfi Prao'V ^Sitfi^ffl fai, 
7-r-K'<**«»f 1. hi. gr ltc<fc0, JWT0D(3 
3)^tcfi£-o r * n * h ;Wliia2&a8feB15maS5EfiiQ a 1 

*#fflrs. fir, »tBL/c^n^ h)\smm&mM& 
tnmwtmm q a 1 ^xr<o o^^comm ( o 

g(Qa)maxknF<£>©ffl) CC^fa^i^TXP 



0 ^Qa ^(Qa)max 
[007 1 ] C<D*&tCtel>T*>, IMZMJzEGR^Jl 

p m ifcm^ *x v ? y'mmmc <t o l fcffi* jbi > r 

*>H^ e K^EPmi l^JfflJ^micD 30 

* F;l/7^^aX-^(iSa f . aJ <hJcJ:otf 
*Xn^ h;Uiia^aSjEa^t*CAQa)inax*KS 
U , C (DS** a * h JWI»S»8SKf bS (A Q a )ma 

ilil^uSS(Q a )max£ffli,>T felt^ 
(Qa)max=Qa c .n + (AQa)max •••(35). 
[0 0 7 2 ] fcLbKJ;^ _hlE ^-Tv '?'S20QX<D F / B 

ETC}g^gP3 7C7)MSi It, ±IB*-r *:?"S200T|£ 40 

x-^lSa^gfflt^ SfctC; *f*?:/S220 
t. EGR*B^»3 8<D*aatUT. ±IB*^:/S200 
•CifHl/cEGR^l/yia^Sl^glQe iv- 
^l/K^EEPmiCcS^tiT, v > ?"#B§GC iOEGR 
^l^gSsffiS e*SttJU Xf-^^S230^it^ e 
[0 0 7 3 ] * f * V'S230X\^ Iff 2 (Onm • *K*£*Uffil 
3 2 aO^IMi LrS^W^S**^4««S 
Gf'^fl, A-f>Srtt»5. C©*»«ttS* 50 



Pmo)-(l - f a i ) • Q e 

•♦(33) 
-(34) 

JK£M»SSf»Gf J£lTOC36)5SK:^-rJ:5Cc, 

<Dpmmti$:Mm vximm<oT ? ^ ^mt^mm u tc 

J£§1*£ rfiJ±*T 4 C £ rt*C * S . 
Gf " = Gf 1 -(36) 

^l^^ff^fflP m o % mmttt&3&& f a 
i\ RSlKfiMKd. B*fflBltABFT*fflC^ 

Gf = d- Pmo-fa iVABFT -(37) 
[0 07 5] «±©®H«Hl^-^>6cStL % 06©^ 

mfgfWSfJSP2 1 OCiSKffliiLr, ^^>^ft-fe>1f2 

r^ttoaRfS^WiJ-r-BAMISfft^ XT*^^*s3icn?. 
^^>^fta«SSP2 2«:j:*^^>^ftS4«JiJ«ffll* 

[0 0 7 6 ] B8(C5*T<fc5tC. *Jf$S|-C«, S^tSO 



(10) 
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BTDC97" . 65° . 10° C A4?&t ^7 > *ft-fe 

#3MlOBTDC97' iSJO^^fHr* 
6 # 1 StScDB T D C 1 0 ' 4 OPe^r 3 l(D«flffi|^ 
#4l@f©BTDC97° 4«r©*^C*l«r** 
#2^I(DBTDC 1 0° <fc©Par2fflcr>aiBf*!lglJ^VI/ 
X, #1, #2^i©BTDC97 < 4fu<D<^fUgJCD 
BT DC 1 0" iOlfflr 1 fflO«lBPPIS0^VU^*sa*S 

[0 0 7 7 ] fi£oT, aiB«Bi|^^^A*Sti-5ff 10 

3WA*Stifc«0^^>^^^bXtt, # 3 xVM<D B T D 
C9 7" ©^V^-rOlO*, 2ffl©««*iJSip^^}&SA ' 

• lffl(D»1S#!IJ5!I^^U^3&SA*Sti 

l«tSf*4C*tt# 2m©OB 

TDC9 7- ^7>^^xt$,i m<D$mmw# 
titf#2««4«w-r*. 

[0 07 8] ^/c, BTDC65' . BTDC10' CD 20 
^7>^{4gte, B T DC 9 7 ° ©*5>*rOUX*6 
©^l/*«T«j£U # imt§fOBTDC9 7" (0^5 
>^;l/X*0 (»ieffi«> 4 UT £^>*-rOU**JA 
^SnSSCCl. 2. 3. '"ii**9>h7^U- 

[0 07 9 ] »<^7 L v^ r S320"C«, #5>^ftfi^U 

@<D^^>^^XA**-C<DigiaB$ra^ft^U. BT 
DC 10' ©*^>^\rt/XA##>6BTDC9 7' CD 
^^>^^*;l/XA^^-CCD^^>^HS9 3° ft<Dm& 
B$fffl*MT9 3 % BTDC97' CD^^>^^VUXA^7 
^6BTDC65 e (D^>^VUXA#£T'<D^> 
*ftl£3 2" #©ffiM*HI*MT3 2 t BTDC65" 
^7>^A;MA*^f)BTDC10 8 ®*-7>*,re 
"^XA**r©^^>^fttt5 5V»©aai^KI*MT 
55iir^y(c^h7T5 D §SiiMMT9 3. 
MT, 3 2, MT5 5<D£fWl 8 0" CACDMiS^ra 40 
4 Ltx>2;>iR«N e (onmrn^htiz* 

[0 08 0] Xf-; ^S330T«, «I*B$8BSJ£S|$ 4 1 % 
*^C««BKS»4 3©«H14m^ «*fB$88* j&U*Jtl 

«*fB$»!T i n j ?&£Afc4^<D*7:>?A*>6 

^-^>-CB3estifcja*WBT i fir*. ^aHEAfctt 

[0 08 1] tit, Xr^^S340"C, 
SB4 2CDMJ14 Lt v ^(DZ^^ZftmLfr&^ltoM 50 
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^ 5 > * ft SB£##^a«£»fc»5£0D 9 5 > i'flgtcfc 
tt£SBi#-e*>S4* % «£A>^x#£££^*-fe* h 

ttrOl/**W>5>**dr 1 0K:ffi^snrjlRf4*si«ftfS 

ft, ±13*7^ ^"S33(rc*s u ft: £ ^ > #-e&jk 

[0 0 8 2 3 W±tC<fcf3. aax>$y>h;b^{c»JCBb 

^*o«a«aw*^s»r*iR©»AS«iRcFE gr 

Xd7 h;U*<D»ft*S*4 0r£D£fQg<D!RA££UI 
&c*fT&:7 — K^c* ^%?f5C4-C. h;U/^U 

^CDfiffiSlftfc: J: £P§PffiScD^ftfc: J: SftKKfL 

[oo83] m&m<Di&mcfcc, mmm 

Mlft^tSiL «A^fiS<D*®{fc, EGRSCDfi 

^^3>(D{£®^^C4^r^^44fe(C, Xh-f* 
*«««cR 'J - >Mi^ fe#»Tj2:fitt2ai!8J:b-c<o*(l 

[0 084] HOawai 0W*«W<^|Jfc©»2JB« 
5C«t>0 , H9 ttjjRfi- • »Si • E GR8«gB<D:/P v 9 

m, m i o - • e GRpmmmfr-rxDy 

[0 08 5] *JKS8« % mm<om 1 fltc*fb, 0 9 tc 
^•Tcfc^CC, »«• EGR»«BI3I53 0CC*j^r % 

»l»«OF/B*WBII36, *2<Dftfir-j|)8«WWv 

3 2acDiniirt^^ I f^jEL, -en^n, F 

2b4*r^44fctc, »UI5«©««lv*-^+ 344 
»l©»«'7*^t34ilyt, COSKD^Ji 
7*-i/t3 4«: l IS2CD»aSllS8I-7*-^^3 4 a* 

[008 6 ] -Tftto^, **«rtt, ^2CD^$fJta)v 
*-S?ir 3 4a KfcOT, #ffi»SR^»ffi<Dffl»W«cE 
*lSS^»Jffir AS^^flc^EE^aWffi P m e e \ 
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s^^EE©«»«ttE*«s^ffl-c*ss»aE» * 2oM • mmm**-i?+ 32a icxzmmfrbn 

M?I«Pmo^fHiU, F/BJW»S3 6 At, 2 ©ft*? * j^ffHjffiJv*-^* 3 2 b (C<fc5«BIK:£M 

EGR/^b^aa^SKMlTJSBH^ffiQe^ffUi-rS I/, * :/S22(XDE GR}gS»3 8 tc<t & E GRM;U 

HL #2S«jaE»»EE*aa« Pmee'i #2*S##EE ^Jg^ffifftiJ&Si Xf-; ^S230©» 2 <7)Jt£3r • jtBjftM 

SaJSEfl^BE^Wffl P m o * t £im&ftE.imm P m WkVt*>*Tv y S22 5*»A"T £ o 

oiOKSlOWFlH8l»tt*ffl^4J:5CCLTC»S. 36 [008 8 ] ^^t? ^S200©F/B«IIBI»3 6 AW:<fcS 

^EfiliaPmo^ffi^tMEBftl*^*!* 10 Pmo, #ffiAti^EJtSil P m e e^fStfflU 
Gf'fcJWBTS. ft EGR/^^aa^SttftJJtBHjetfQe 1 . XD 

[0 0 8 7] *JB»rtt, CROMER • ^ h^Wiil««aHi»3JBRjefaQ a 1 ^gurrs**. C 

KM - EGR©Jt0M3;t^^>&t#U SI OfCSTJ: <DP£, *lJBJBCC*tU JWTOC38) , C39)^CC^*T «fc 5 

*>&c, ^^»^s2<»ecfews«a©rt8*F/B»iiw» ec % -en-en. 1 nmmmm<o^^iSu^Emm<on 

3 6CC<fc5^3^6F/B«IIBSB3 6ACCJ:S«iffltcgg r^ffl^HI I me e ( _ 1)s 1 MffJHniiroSftJSfi^SM 
JE-T & £ t 4> CC . * ? v ^S230tC*J^ £&SOF*3^£ * ^©B#Pb1«^MI Imo t . n »6 e 
Qe 1 = h2 • Pm e e + f 2 -(Pmee^-Pmee)* 

g 2 • I m e e c-i> --(38) 
Qa' = hl -Pmo+f 1 • (Pmo* 1 -Pmo)- 

(1-f ai)-Qe + gl • Imo c . 15 --(39) 
[0 08 9 ] ±12(38) t (39)3;-C3?m Lit EGRa';^ ^EllffiPmo, K;Uffl«^8fc»»SffiQ a % 

•aBft^M£«J7BHKJEfflQe\ * n * K Jl/»»SSt«£ EGRM;u^»«#x^t^tf[Q e , ^fihtJtSfflf 

■KzjwssfiQ a 1 mi mmtmm. o ^e>*^ss ai, 1 mmmmm<o?^m^smm<omrAm^m. i 

SO©fflfc:g§fD3i*\ EGR^^^iiia^SHiaSffi moc,,. 7 - KM ?tf$k f 1 . hi. glfcj: 
[0 0 9 0] — 2f % *^*:/S225te:foW*»2<D»»«l Bafflr*4^»S»»ffBaffijEffiPmoh**JMT 

Pmoh'=(l/f 1) - (Qa + (l-f a i)«Qe + 

(f l-hl)-Pmo-gl • I mocn) -(40) 
[0 09 1 ] S *^U5##EEfttSffiPm e e % 30*R^^U^aift**^^«(CtH^TS^aaffl"C*S*ffi 
E GRAO^IM^iSIiSiQ e . 1 ©J®i§R8m)(D SU&^EEaa?§IEtfIPm e e h'*JWT©C4l)5£K: <fc 
*SateR»»BEK«©«FlH8l»fll I m e e c . 15 , v j - -.tfffiT^o 
KM**flHRf2. h2, g2K£D, S:7E$n/cEG* 

Pmeeh , = (l/f2)-(Qe + (f2-h2) - Pme e - 

?2 • I me e c -i)) -(41) 
[0 09 2] S«Jlfe»»BEB«l»SjEfflP-m oh -fc (42)S«C J: 0 % Zmj&frfrE^mmP m o ' *»ffl"r 

Pm o * = ( 1 - f 1 • ba- d t ) * PmoV,,+ 

(f 1 • ba- d t ) • Pmo h* --(42) 
[0 09 3] #£«US##EEB*»IEfBPme e ♦flMUfiffttbeftffiC^ «T©C43)S:k: J: 0 . 4¥^Mj& 

h \ 1 «t»H JHWr(D#KS«fiS»»ff : F»Jffl P m e e * OAE&Wm Pmee^IWn, 

c-i,. a«Jt*SE«f a i. y4-F*v*mkt2 % ♦ 

Pmee'=(l-f 2« fai • be* d t) - Pmee , ( . 1) + 

(f2 • f ai • be-dtJ-Pmeeh' -(43) 
[0 0 9 4] fit, ±E(42)Sr#UOfc£«JlR## * Pme e'i^»MiSIPme e itJORi© 
ffi^«fflPmo - ta^»»E^«Pmoi©Bai BHRfflt&fB I me e*, JaT©(45)5£fc: <£ otSffit 

ii^C 4 ±IB (43)Si'C»tti L d^BE^Jtf * 

Imo =lmo c . 1) + (Pmo , -Pmo)'dt --(44) 
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Imee=Imee c _ 1) + (Pmee*-Pmee) - dt -(45) 

[0095] ffiswccra* jjbc40)5«(c isffiaj^^ uta-r&o 

EB«tSiEfiIPmo h\ ±f2(4i)5£*c J:S#ffi«sK» Gf = d • Pmo" • f a iVABFT -(50) 

ftmmmmiEmPm ee^it ^-n^n, j£nF©c46), [oiooi c<ds^k^i*»«hig fcoga 

mo'\ #3gafiR^BEB«ffi«nW»5E<iPme e M i ^t3 2 b &C<fc£^]&5M*JI^$MBPm o 'Zm^T 

Pmoh" =Pmo H -(46) 5^C, ^fl^^BE^aiffiPmo^fflli-C^OSlJfBl 

PmeelT = Pmee M —(47) io tSffaccW^S^ffEE^OjCSfiBSattWCC^aj-ra 

[0 09 6] C©»^ v 6 ACCj:£jzyi Ci^Ct^ X n ? h;^-^E GRSO^- Fttftftft 

cc*$i>r . e GR^^^aai^iun7«i«3e«Q sn^stw^noiincc j: ^rnB8o8Mc0Rcc^6 

e\ xp ? h^mm^^mmm^mQai^. ^n zmti&mmi,. mmcomwrn^m^ ox^^Mcomm 

-en* o^6*^ci»«©fiH«c«a«i3-i*SR, *zrai$5E ^^ttc*f*rsiisett*roj±L, ^jBSJtffjfecDiwfiiK 

e, A^r-fe^ h/*yr-r*<fc5K:u »2<D»asoi» [oion g/c, U2©mj&scc*f 

'7*-Pt3 4a«:J:4«I'C t «l^7 ^ yoffifctS U xn v hJV&com'&femt LX&C2>mm<D&Xg. 

DT2«flt»»EBRSI©«FlHa»ffl I m o . *ggSUfc# SUKCjPtT*:? -f - F'<* ^£*t 5BL EGR/^I 

#EEsRM<DB*lfflS#«l 1 m e e ^SS^T-S J: 5 Cc Lt *8fc«*^ISffi4IQ e 1 SWd * h^SEftffitujfE 

M< . WJ®fBS**S^5B«UooltS©^f*««-r 20 IfiJMSIQa 1 ^, WgiCDS^ffil me e, I 

[0097] -r^^>^, EGR^^^iiia^/xgfcat7]«i ft^frjhu F/B$fjwss^r^±*r^c<h^T^^ P 

l^ffiQ e 1 £ E GR^^l/^ffiift^SSaS^fitQ e <h [0102)01 1 S 1 3 3 

^uc^arfifD^^^E^^yrL. egr^i 5c«*>o. 01 i«^>^>$iJtsp^©^(t«^0, si 

^^^^^^0 6* <hEGR/^;U^ajg^^ 2 tt«B4 • • E GRfrjffllgpcD^n ^ H 1 3 « 

h^mM^^^&^MBQat^mL^t^mu SSfcStU SI lK^frTcfc^CC, y&gVgffi;fr-k>1t 5 

iQa'iXP? h^3lfi££^fiSSS^fiiQa£^S^ 30 &££&K:, X P ? h;loliaS»ife»*9f*8BTS»A^ 

[0 09 8] ^LT, fef07 7^A, E*5*CC^UTS 2 0*Cfalr>T, i»a«EE*-fe>1f 5©ffl*CCS-^l*r ^ 

hri»Si*, ±8B(44)j£teJ:oTS^##ffR^cD -*^F^ffiPm^»Hi-r^^^^l/F^ff#tBS|52 6 

B*IH«#flI mo*»lljU fflft^^A. EOC^fh 3\ «ffiH2>*9<Offi*(C»^C»r^*^ F^ffPm 

fr-^-fetrhShtMitKB, J!aT©C48)3C{c^ *ffW"r*^-*^ F±£E#ffl»2 6 AficSHly, R£i 

t^^fc, SMA^SHieo^lini^a I m o * 1 m t«cM3hfc!RAaJttt>f8*6©Hl*K:I^» 

B^H^ffi I me e*#tHU fi&fP^^^E^-fe * F § [0104] g/c, 0 1 2 tCTjVT J: 5 &c % • • 

nri>£<fcSNc«\ £iToc49)3SK:^-r<fc5cc % 40 egr*jsp8B3 otc*jc»r«, F^aia^asaat 

<^£- (is, Hi2tt % m2<owmMm-?*~t>+3 

Imo =Imo MJ -(48) 4a*dto»2Jg«CC»L. F/BWWBP3 6A*F/ 

I me e = I me e M1 -(49) B $00)3:3 6 B^Ml/cM^Lt J&20D5R 

[0 09 9] ^^^^S23(XD»2©ftW-««IW t 3 4 a^L^cl^ 1 ®i«:Jtl 4 F 

«pv*-y+3 2 bCCcfc*«fflirtt, Jg2cD^»^ /B*»J»W3 6*f/bwj«s»3 6BccaHi/r*>a^ 

+ 3 4 a-CffaUfcgajS^EE^ajfflPm C<h«*7^r^^ 0 

o " , SMItRSffl f a i \ R^KffiK d . SffiSjKJt [0105] *^3Brtt % H 1 3 K^iT cfc 5 CC, -^B# 
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comic. SS-b>it9©ffi^icS^t^rv^^;l/ F±EE 
JitfjCCg^TXn ? h^ilil^a&MttS'JfilQa v e 
[0 106] tftto6» x-r^ys7iro-7^^;l/ F£ 

«T<D(51);^C<£^ ^^^^K^JEPm^rgL 
fflU */c, Xfy ^S72rc0XP y hJl4U££ttHI * 
Pfo = (l-a-dt)-Pf Of. 

(ba- 

[0 10 8] *©ffiOKMtt, »l»«**l>tt»r2» 
[0 1 0 9 3 

*ffi«flE» & *«B£ 0 6 »««WSCf EGRSM 
fflr^ lo, xn* 

©»A2S«*K:SFt-r47-f-K^^^^T ^ttf>. XP 

©^-FWttJ6SHBn*»flir*i<!:«>CC, XP*F;l/ 
h;U/^uy(7)f@^CD^^^ + CCJ:a?SS<SM^. xp 

7HtMt*sfp&fts. 

• EGRSlJffllSBcD^n **0 
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A/D£»U h;uajftffiflBEJBf«fflQave 

Pm = p • R a - Tm — (51) 

[0107] fit, »l»JBOF/B»l»»3 6*4 
C 2 J&m<D F / B $»J®SB 3 6 A CC *5 W 2> G6)5ttC <t 
SffilWK^BE^^mP f o©JHH«:a*U 
F/B*llfflJg&3 6 BCCiSJffiffl-Ctt. «TOC52)SCtC7r: 
*T «fc 5 K . 1 ©JSPSSBmrcDX p ? h ^iliB^MSKSQ a 
io c-,j«:«^r®A2»«H2>if 8(Dffi*ccS^< xn * 



_ 13 + (ba* dt)-Qave + 
•dt)-(l-fai)- Qe ( . n -(52) 

[S3] [*Lb «3«!ib^-^>©7n-^ + -h 
[04][iJ_L SS9Ml^^>©7q-^t-h 

[06][eJ±, *^>*ft«0ii#JU-^><D:7n-^ + 
20 — h 

[0 8] |eJ_L «mflJSU©lftBJi0 

[09 ] *mi<Dmm<Dm2Bmt l c&t>*) . mn • ®m 
[010] r±, EGRmmsm^-^> 

[011] 2^«IBa«B»<D*3^JBCc«*>D. x>i>> 
^Ji»^<7)^^^0 

[012] RLb. 8**4 • «a • E GR^fJiaigpcD^n ? ^ 
30 0 

[013] |hJ±, SW^^-^O^a-? 1 *- h 

1 -x>>>> 
1 b-xn v Y)\,j^)\,rf 
1 >^*** 
1 3-XP ? F;l/T^3^x-£ 
1 

Te-B«x>^> F;l^ 
Gf ---S^ttftigttjt 
40 f a i ~i!ijWt*Sffl 

Pmo •••ffijJW8^BEW£« 
Pmee -i^Msm^Wmm 
Q e ••• E G R^U^jlj&tf^S 
Q a -X n 9 h ^iiiS^fufcS 
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